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SUMMARY 


This leport contains a listing and instructjons for Progiam NONCON 
for which the mathematical development is presented in reference 1 Tne 
computer piogram determines tne distribution of lift, pitching moment, and 
pressure on nonconical slender bodies with leading-edge separation Input 
data is obtained from Program SMITH (ref 2), which is valid for conical 
wing -bo dies only. 

The program is written in FORTRAN IV language for the CDC 3 600 
digital computer at the University of California, San Diego With only minoi 
modifications tne program has been adapted for the IBM series 7040 and 
7 090 computers at the NASA Ames Research Center 



PROGRAM DESCRIPTION 


Pioj^iam NONCON cons ists of one main pjot’iam conl.nninp 
lliree subroutints 'I he function of c<ic!i subioutitio Iia<^ been ‘^piei(i< <1 m the 
LibUng 


Input Description 

Piogram NONCON is designed to read two kinds of input data 
1) regular input data at the initial station from Piogram SMITH, and 
Z) restart input data for continuing the piogram from any station at which it 
may have stopped prematurely due to the time limitation 

The initial input data are read in by using the format-free subioutine 
INLIST provided by the UCSD Computer Center At most other compute i 
centers and with only slight changes, the data can be read by using the 
format-free NAMELIST feature of FORTRAN IV version 13 


Regular input data specifications . -Definitions of the various para- 
meters used in the regular input are given below 


Configuration and sheet parameters: The following data are read m 

from the Program SMITH output. Units and normalization are the same as 
for Program SMITH, viz, , velocities by the axial free stream compon- 
ent U, and vortex strengths by U tan 6. 


NS H APE 
RADIUS 
A Z AXIS 
BY AXIS 
N 

BETA 

DELTA(I) 

AE 


for circular, NSHAPE = 0, for ellipse, NSH APE = 1 
radius of circular body 
semiminor axis of elliptical body 
semimajor axis of elliptical body 

number of sheet segments, number of pivotal points 
(3 < N ^ 25) 

dihedral angle 

difference of polar angle between two pivotal points 

desired value of angle of attack parameter, 

AE = sin a / tan 6 
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ystarv 


real cooidjnale of isolated vortex in tiansfoimcd piano 


ZSTARV 

GVORT 

D(I) 


GS(I) 


imaginary coordinate of isolated vortex in transformed 
plane 

rsolated vortex strength, normalized by U tan 5 

distance from the isolated vortex to a sheet point 
I = even for a pivotal point and I = odd for a middle poii 
only the pivotal p^oint polar distance must be read in 

sheet strength at a pivotal point (I = even) and middle 
point (I = odd), normalized by U tan 8, only the pivotal 
point strength must be read in 


Nonconical shape parameters The following parameters arc defined 
in figure 1 for* two configurations, either with or without incidence The in- 
put of Program NONCON must be'inodified for different configurations. In 
order to read in tabulated data for nonconical shapes, the main program and 
appropriate subroutines must be changed. Any set of consistent units can be 
used for the following input variables, except that angular input is to be in de- 
grees. The axial component of the free stream velocity U is taken as unity 
in Program NONCON, and all output vortex strengths are therefore m effect 
normalized by U. 


ICASE 


ICASE = 1, double-delta wmg configuration 
ICASE = 2, ogive nose configuration 

LING 


LINC = 0 without incidence effect 
LINC = 1 with incidence effect 

XO 


initial station and starting point of nonconical section 

XF 


final station of the configuratipn 

DX 


increment in x 

XCF 


final station of nonconical section 

XCFS 


station shown in figure 1 

spanf 


semispan at the final station of nonconical section, it is 
sufficient to know either XCFS or SPANF 

DELTA 


semi-apex angle of initial conical section {deg) 

DELTAF 


semi-apex angle of final conical section (deg) 

EDETAF 


final incidence angle (deg) 
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Tolerance limits 


ACL upper tolerance limit for all three loops, if the oi ror 

exceeds this limit, the program will stop automatically 

ACCl limit of percenCchange in vortex strength (occurs in 

LOOPl) 

ACC2 limit on adjustment of zero force condition (occurs in 

LOOP2) 

ACC3 limit on percent adjustment of sheet shape (occurs in 

LOOPS) 

Iteration limits 

NLOOPi limit for LOOPl, generally set to 25 

NLOOP2 limit for LOOP2, generally set to 10 

NLOOP3 limit for LOOPS, generally set to 15 

Parameters involving pressure- 


LPRES 

LPRES = 0, no pressure will be calculated, LPRES - 1, 
pressure will be calculated at prescribed station 

NP 

total number of pressure stations 

NPRE 

number of locations on the surface per station at which 
pressure coefficients are to be calculated (NPRE< 100) 

PDX 

increment in x when pressure coefficients are to be 
calculated. DX/ PDX must be an inlegei- greater than, 
or equal to, one 

PRPRIT(I) 

X station at which pressure coefficients are to be calcu- 
lated (I^ 10) 

Restail 

parameters 

IRESTA 

IRESTA = 0, program starts at initial station, 
IRESTA = i, program starts at the station wheie it 
stopped previously 

LPUNCH 

LPUNCH = 0, no punched output, LPUNCH = 1, punihed 
output if the program stops due to the time limitation 

TIME 

maximum running time 
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END, STOP caids An END caid should be msertcd aftei the in]3ut 
data fox each scpaiatc case After all cases have been specified, a STOP 
caxd IS needed to end tfie reading piocess. 

Typical 1 egular input data* Typical regular input data aie shown lu 
figui G 2 The Older of the data can iJe alternated. Some of the unnecessaiy 
data can be neglected as shown in the second case of figure 2 


Re-start input data specification . - The restai-t input data spen fica'i ion 
IS divided into two pa its. The data Tor tlie first part consists of the same 
input parameters requiied for the regular input data specification. The serom 
part consists of a set of prepunched data cards (with a definite formal and 
order) obtained from the output of the previous intcirupted calculation. 7'his 
part must be added after the END card of the regular data input as shown in 
figure 3. 


Output Description 

The calculated results will be printed out by the program at the initial 
station as shown in figure 4. Print out at intermediate stations is similar, 
except for omission of the title and some parameters A description of the 
important output parameters is given below. 

General output specification . - General output parameters are defined 
below. Quantities referring to geometric properties of the model or to 
properties of the vortex sheet in either the physical cross -flow plane 
(Z = y -t- 1 z) or' transformed cross -flow plane (Z* = y* + i z ’' ) are defined in 
figures 1 through 3 of references 1 and 2. 


X 

S 

A Z AXIS 
BY AXIS 

YO 

ZO 

GVORT 

YSTARV 


X station of the configuration 
semi span 

semiminor axis of elliptical body at initial station 

semimajor axis of elliptical body or radius of circular 
body at initial station 

y coordinate of the tip of the wing in physical plane 
z coordinate of the tip of the wing in physical plane 
r^, vortex strength of the isolated vortex (based on U = 1) 

V , real coordinate of the isolated vortex in tiansforincd 

V 

plane 
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ZSTARV 

RSV 

THETAS V 

YV 

ZV 

RV 

THETAV 

H 

I 

GS(I) 

YS(I) 

ZS(I) 

RS{I) 

THETAS (I) 
D(]) 

Y(I) 


z^, imaginary cooidinate of the isolated vortex in tians- 
formed plane 




1/2 


0* = tan“^ (z'Vy") 


V V 


y , real coordinate of the isolated vortex in physical plane 


z , imaginary coordinate of the isolated vortex in physical 

plane 

/ 2 2 \ 
r =(y + z ) 


2 2 \ 

' + z ) 

V V V ^ 

e_^ = tan ^ (z^/y^) (deg) 


polar angle of each middle and pivotal point (deg) 

index for vortex sheet points (I = even for a pivotal point 
and I = odd for a middle point) 


strength of the sheet segments (based on U =- 1) 

iJU 

yT, real coordinates of the sheet segments in the trans- 
formed plane 


z"^, imaginary coordinates of the sheet segments in the 


transformed plane 

r* = fy*^ + z*^J^ 
1 ' 1 1 / 


0*=tan ^{z’^y*) (deg) 

polar distances of the slieet segments 


y^, real coordinates of the sheet segments in the physical 
plane 
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Z(I) 


R(i) 


z j imaginary coordinates of the sheet segm'ents in tl.c 
physical plane 



THETA(I) 



0 = tan 


(z^/y^) (deg) 


Additional output are punted in self-explanatory titles, such as 
ANGLE OF ATTACK, LINEAR LIFT COEFF, etc. 


Pressure output specification. - At prescribed axial stations 
PRPRIT, upper and lower surface pressure coeffients are printed out 
Pressure output parameters are described below. 

Y real coordinate of a point on the body or wing surface 

Z imaginary coordinate of a point on the body or wing sur- 

face 

CP pressure coefficient 

* * = 0 reliable data, =!= = 1 unreliable data due to the un- 

defined log 0 


Additional output - If the tolerance limit has been relaxed inside a 
loop, the maximum tolerance will be printed out before the general output 
as shown in figure 5. 

When the indicator LPUNCH has been set equal to 1, a set of data 
cards with prescribed format will be punched out if the progiam has stopped 
due to time limits. This set of data must be kept in order, and will be used 
as a part of the restart input data 
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PROGRAM LISTING 


The FORTRAN IV listing of Piiogram NONCON is given on pages 
15 to 53 


Air Vehicle Corporation 
San Diego, California 
October 1969 
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XO XCFS 


I) Double-delfa wing configuration 



2) Ogive nose configuration 


Figure 1. - Nonconical shape parameters for 
double-delta wing and ogive nose configurations 



azaxis-''. ryaxig=o. n =6 nriA-n, 

'■ir L T A|(= . 1 '-.'W' f ^t‘ = . 0^1 <) 1 'lOo/ 7f 1- P ! 

T ^A't ZSTAI?V- . T'f)! BAA 1 M fA/OpT = i ? . 3 1 7 407 

'.?] AAA' p/< ,c. .l,'Al'’Z'‘’ippBP ,0 . 4] 98 AA71 

fM?)-.; ,3''' I ' ♦«> . 34 'jBi"W3')16 .0 , 1 77 

V '' = « ■ 1. 1 lAfiAAAYfx / ) ,0 ,«'?7^A‘S' f,-,l 7 ,0 .61 90983 Oa'-*'' 

:j"(n = «‘' . A ] 7 '9^p / ^ 3 ,.i , A '^78 J A j /t ^6 .0 ,4878 37' ''7i-,A 

IC.A6^=3 I ni'" = x Xv.-op.T xr = 4fu DX = .7 xcr = 30 , XCFa = ? 6./, t;nA„c-- 
riFLTA = 6, nrLT3'='^id. fdfltf=^. 

Aa=].r -7 ACC = 1 .C-6 ACC? = 1,C~4 ACC8 = l,(7-4 NL(CPl = 9p^ r’l (^9l>7 = i {> 

L9pF-- = ] hPRF = 6''). PRPPIT=19.7 74, 7P.B 3?. 6 38.4 ,MD = s P 0< = .7 

I .FSTA^"-! LPHNCh=l Tlf'-t-B. 'M-THOD^l 

Forcc=o. 0. r, , 0 * o. 0. 0. 0, 0. 0. 0. n. 0, 

r.RD 

L[NC=0 ICASE=? X0=3.4 XCF5^-14.4 XCF=12. XF^33.6 nx=l, 

DFLTA=1B. 0FLTAF=6. 

A c- , 8 34484 1 3 76 33 1 

YSTA9\/=c 1 6797^70474 ZSTA^ v= . 6 171 4 0 1 1761 GVORT^l .47501 OYoo-^io 
n = ,396055'= -^799 .0 . 26797428-=90 ,0 .70531561 1 14 

n{7l=.0 .17188580743 . 1626531711 8 .0 .13860561581 

036=. 0 ,71351847794 ,Q .17973149339 .0 .11870438867 

GS(7)=,o .08873588095 .0 .07641131677 .9 .06985709791 

LPCFS=0 IREx3TA=^^ 1 IME=3. F-FTHO'a=2 

rORCr=,265 .6199 5*836 .5183 .402 .3^367 ,794 .27? ,758 

^F0RCF( 10) =.451 .398 ,819 .770 .754 ,751 

’*^MD 
STOP 



Figure 2. - Regular input data 



r'li'lip; P/'l (i(> AZ^X]i!i“0* RYAXIs^ — 0* ^i“6) BLTA"0» 

f'l U Arir=At. TA. iAsAf*'. .f>cAo< Au-,>.''ixiiino^'rr''>? 

vr, ]• fVi?V = . if!.| 1,-'1 A . 1 /.STM’\ - . /QM 19 GVORT^l ? . 3 19 7 ?9 ] 7 /+q 7 

n--,f' . 6 A 'I? I A', 'I ",19. .0 .'-Or, 1 ^,= 9U?BG .0 .A198BA7]9?R 

h(?)-.f ,3 7 i‘>< jn;,r .37iO? 18 79916 .9 ,13,19860917/’ 

] , 11 Y9(wr) I'"-' ^'O -,0 . 137''*f-.9P7,61 ’ ,0 .61 9099307 ■’>5 

C'MVl'.f' .61.-^' 967^U63 . 9>3 2 B 1 9 1 A 6 .0 . 9 828S 7 9 2 26 6 

ICA69SI |IIIC='' X 6 = 22.J Xr-AH. DX=.? XCr=30. XCi-S-26.A 

ACL-'l.r-2 Arc ]=-] .1 -B ACC2=l.L--9 ACC3=l.t;-9 i\LOCPl = 2'J> NLOOi>2^ir MLCrP’'--’ 

OfLlA^O. Uf-.tTAr = lB. P.DLLir = 0. UPriE5 = 0 Jpr61A=l LPUNCH = 1 TIM,':/'-. V.trP'^'r' 

r.xD 


STOP 


2 2 
3. 1662679938 ] £ 00 
2.40A39667164E 00 
?.n?4383 00633F-0.1 
2.1421 874S92HF-01 
1.3B995049413F-01 
1. 239712795978-01 
2.50664184U2 7E 00 
2 .885350780628 00 
3.1032630O270E 00 
3, 109SS992294F 00 
1 .02194302756E-02 
1.T0160485365E-01 
3.41317538783E-01 
5.04786108865E-01 
1 .88632948560E 00 
1.23181086106E 00 
9.63525836372E-01 
8.44755582907E-01 
1 .46866133437E-01 
1.3400871 35-76E-02 
5.65452754009E-03 
5.21276133077E-02 
3.33799156942F-02 
2.32000000001E 01 
-0 

2.46110541880E 00 
3.0498372291 8E 00 
7.5963772729TE-01 
2.1?282n68114E-01 
5'.20151248854L'-01 
8.25902452855E-01 
7.55?6fl65504?E-01 
2.1189817899 i^E-01 
6.5541533298 18-01 
9.7729174851 OE-01 
7.55268655042E-01 
2.11 S98178997E-01 
3.45106935920E-02 
1.13367219199E-02 
5.2 1446796786F-01 
8.2 1586544398E-01 
2.08258352656i:-01 
3.3061907069 lE-01 
2.440860652131 00 

I.2201507J063E 02 
0 
0 


Figure 3 


9.36755289A65& -01 
2. 1550CT81371 IE 00 
2.85843133437E-01 
] .78 2229934 5 OF-01 
l,3144oon4537E~oi 
1 .24Q78599535E-01 
2.63388428564E 00 
2.985^80fj25Q3E 00 
3.12379829149E 00 
3.06606229511F. 00 
5. 17073093844E-02 
2.26430768432E-0I 
3.97226297995E-01 
5.572S5749062E-01 
1.61702954397E 00 
1.12679983323F 00 
9.054628673niE-01 
8.27621857356E-01 
4.0379K)73C35E-0? 
7.74238613294E-03 
5.859131651 93E-03 
4.38515639689E-02 
3.30795402988E-02 
2.00000000000E-01 
-0 

2.72830461484E OO 
3.078839n3506F 00 
5.83741734060F-01 
-4.22586470225E-02 
6.24075876630C-01 
8.353892n331E-01 
5.84120523708E-G1 
~4.20032555?26E**02 
8.11666935240E--ni 
9.9911747383'3E-01 
5.84120523708F-01 
-4.20032555226E-0? 

2.85879631684E-02 
-2.27140075742F-03 
6.25742809975E-01 
a.36970259496C~ni 
2.49890707544E-01 
3.34303377503E-01 
2.82860181 550E 01 
1 .2435481 2180E 01 
n 

2.4408606'- 21 3f. 00 
7. 00694 90 J 152F-02 

Restart input data 


1 . 94 14305 099 6 E 00 
1.06678235892F, 0 ^ 
2 .5n21449f>409F-01 
1.586005:5143'-,F-0] 
1.238639‘^0659E-01 
l.?Al85Q1947ir-Ql 
2 . i’7725470049E 00 

3. f'487oa59755E 00 
3. 1 32688194 16F 00 
2. 9914641 1393E GO 
1.04467877524E-01 
2. 86376671848E-0I 
4.51 1979862706-01 
6.973634060918-01 
1.42442020250E 00 
1.021788805438 00 
8.75109225116E-01 
8.10488131805E-01 
8. 36823538691E-02 
6.39998454973E-03 
6.52106430603E-03 
3. 69191332720E-02 
3.404181 16684E-02 
1. 59758437600C 00 
2.27999999997E 01 
2.93347665429E OO 
3.01312555937E 00 
4. 2248582?8CrE-01 

-3- 503 R4Q34388F-01 
7.490522410608-01 
7. 7665366h692r-ni 

4. ?l56757578‘'E-0l 
-3.405778876 ?6C--01 

9.0679698886PE-01 
9.36907306197F -01 
4.21567575785E-01 
-3. 4O577887676E-01 
2.19045849895E-02 
-1.05256575854F-02 
7. 505918463B1E-01 
7. V78939266S1E-01 
2.99844262726E-01 
3. 1084871 4871E-01 
7.880669027331- 00 

1 . 923'=84 5 ?2 52E 0? 
2.44086065? nr 00 
0 
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. ’ Re al 1 rt input 8 it t 
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NONfoNlCAL SMI'TH 'PROCSav 

. . "circular BODOTrTT's'TRA'tGHf ST rJFes " 

1 

X= 22 80000 S? 2,19638 ""'AZAX'lSi fl BYAX1S= 1.59758 

. - y0'="2,'39636'e "OB ■“ZO’J'1 0“ 

' AWGLR of attacks 18,00000“ degree ■ 

' - - YSTAR’/s 9,?g64AE-01 ““2STARV* 1.9071J6E OD 'RSVs-eiUB^SE TO THFTAS V=“T.'12l22e 00“ ' 

— ■ YV= 2,116^0t 08 Z\f= 1',“049A9F“OIT RVe 2^-36-233F O 7HETAV’ '^.BOSARE-Ol “ ‘ . . - 

Hs~2-,2e5-47fe--0a r-,5-669TF5Tn 6rrff5tI42E-0'I 9TITr3-09rE=TI ^ITfaTT^ E "0 0 rr370&FE“‘DU 

- ' 1,59843b 00 1.82673S-00 2.05513E 80 2 . 2S347E''5'0' 2.5ri82E'0D- “2.74017E 00 

I US --- -Y3- — - 7S rS TKETAS " D T T R THFTA 

-r~l RTArf’-tri — F7l33-="Crr-V,-5-91-ET.-(T2"““37 SSTtEFin: — iTSWFFffI — iTS'5“6TD"i5 — rTTSITTo ?T4S'8'F-'o2 2,461F“tio 

“ 2 " ?.806fc“!rr' 9;320=-Dl” 3;n9E-01 ■9“,834G-(J3.— 3“.249E“01“'T.'5'94E“00"*2:58Be “o"6 ‘ 1 . 32 l 6 '-Br'T. 592 E W5;o99b-0? ' " 

3 2.466t--Ui— -i-,?^R = -00 ■5.<T32E-01 "Y.361E'"(J0' "'47105r-01‘"'rr4r3E"0“0"”2':ri4“E""o‘o""2';83TS-oT"' 2". r2"SE”00 "i.'l)4()fc-0r" 

- 4 - ' 2 ‘. 1 C 5 =-- !rj~ 1 , 460 = TO B“'2i"2e^(ri iT579“E~0'0“ 5TT2'5F-01 lTRi3e~To Tr794T“00 r.T^Se-ol 2,“8T5"e' To'"Tr70'Ob-uT ~ 

5 1 .753“c;i-Ui- 3:.629S"iI0“'“r7O5“4E“ 00 ■ i7"940rTo“ "5;r4TF‘-n‘"'rril'0E’ 00 ' 2‘.8S9E'"o 5" e^SesET.oi" 2.933G 6 0 ' 2 .1->“3EVo 1 

"6* 1 .btge-Ti" T.72 i£“00““1“.'299F 'OO' ■■ 2 ”.'l 56 E“" 0 "(r" 6 V 4 "<f 4 F'F“o“i i“,'n06E‘“’Do 27674'E'W‘ 8'. 4'6'4'£-'ol 2 ' 976E'To" "'27864¥-0'f 

~7 “1 " 3 r 6 d'Di rr 7 T ^="00 rT 5 l 9 F~ff 0 27345e“00 7,049E^6T fr^7E«0l 2 , 874E “otT" TToTIt To 3 . OSbE'To” 5" 4 T 5 b~oT 

-ft i 2-»3 F-ui" 1-.-V95=" ffO "■T,7'36F““60 2,49B¥ Ob'"'?'. 6B7£’F0i 8 .YlSE^W '2 ,~831E' 00 1718BE "ob""3.07"0‘F 00 3.9/5E-bl 

' <5 1 a-ftt-Gi 3.,781'S'OO' l,'/36F"0t) 2 .■ff32"E"0 0' "B , 275E'bT'"6",612E4b'r*'2777TE'’ lfo"'a ,34 4g “oT" 370 r9P' 0 0 4'.5l"«>6-0i 

31 T. 2T"6T •=Tf 17723T OO 2712 ^T'C ff “2 ,T35T1£lO r.TbUTTT 87TneWl"'TT6'^rT6 uVf9i~Qh 3To56't~ Ol” 5 , 052E -ll 

II , r,ft38= DO 2,292E“lO 2 ,Oi7E'00 ' 0 . 5 0 2 E- D 1 ■"a , 13T9E'7&T'"'’2 ,~556E' 0 0 ‘ i7 595E ' 0 C *’ 3 . 0 1 sIE 'O 0 5' 8 /ce-O'f ” 

I _ 

!2 1 "'6t--i l.F2'j= 00 7,473F 00 2.369E'00"1 0i?F" 00 j" 7‘, 96SF-'or“ 2", 413¥”P0”"1 ,6f 9ETo “2794 0l"0 O'" SVn'79E'0T~' 

rtrJFAR-riFT XtrbPfT: -r.39008E“FOr NOWLl'NEAR"X7rr O;OErr*~?,T9672E*0i""'lNTDri;iFf~CltErF£ 0" 

TCTAu Ll-rr COEFFa 3-586S0F-01' ^ - - - 

LiN'EAR ‘10-lF.'|-T=-9,?ft77'£«r2 ^40^i1 INFAR HORFNTs 1746443L«01'“ TOTAL' MOH5M COcFr=" 2“. l9i20E>0'i 

w Figure 4. - Output at the initial station 




-S I'JJ L/.C}'',C 1.1 5745E-04 ACC«5= J . C 15 0 0 OF"0 -T 


the results at X=' 23,20oOO ANU S= ' 2.44086 ARE 

■ * ~ “ ■ gVo'RTh 3;16527E~00 

Y3TARV= 9,56795E>0i 'ZSTAftVs 1.94143E DO' RSV= 2.15563E 00 THETASV= i.l2l?2E OO 

' YVs 2.15501E 00 ZV="'l .'06678E OO "R\/» ' 2 , 4T460 E 0 0' ’THETAV= 4.59659E-01 ' 

2-s RS fTTETAS D~ Y Z ft 'THETS 
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! 
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I 
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Figure 5. ~ General output 



PROGWA'^ NOHCON 


C THIS PROGRAH I, ITISJGNrn TO CAUCULAlt: VORTO SHTCTS BY 0^1111 

c f‘KV'CoracAL thlory 

Cc;H. iOf>!/ 1 /Rau !OS , AZAVIS, bYAXIi>. i. , BETA 

C0M''0,V/?/HSHAPr 

CO: » / I, xxx. X> ■ 

CC)('■^XO^J/i/ALPHA» AL, DLLFA. DELTAFt SPANF* XSPF • lAAAF. SINAI » S I N 
1?, COSA2, rOSAT. COSP.T , 

COMMON/5/f’ I ♦ TWOPl, Prf'2, THPT02* RAD» N« NZ. MPFl- NJZPI 
a-vu/OM/S/OVORT . GVOPTOs GS(SO)* GS0(50), RS(S1)» tHETaS(SO)- YS ( SO 
11, ZS(SP>. D(51)- H(S?), R(50), RtKSG), THETA(SO), Y(S0). Z ( so ) , - 
2SY, THFISV, YSTARV, ZS1 ARV, , TV', ZV , ZVO, R\', RVO, TMETAV 
CtY’'.10N/ //RO, THtTAO, YO,-ZO 
COHMOM/9/EPS 
COH-'AOfi/lO/LinC 

COVMOrj/l 1 /Dlv'DSGR (2S ) , DWUSGI(?S), r(2S), C0SPH1(2S), SrNPHI(?‘>). 
1BS0(2S 1 

COWMQO/12/X, NPRE, IPRIT 
COMMON/] 8/ACL 
COMMOH/20/XCPF, ELELTF 
COMMON/2 1/DEL TAH( 50 ) 

COMMON/22/DGVDX , DYSVDX, DZSVDX, DDOX(25)* DGDX{253» A, XY, Tinr, 
IPPPRIT(IO) 

CO‘'^MON/23/ICASe 

COMMON/24/NLOOPl , ACCl , SUMK25)* SUjV 10(25)» SV.IPUK2S), SSLMl(?>->), 
1 SUM2(?5>» SUM2UI25), oUM21(25)» SSUM2(25)» CSPIX0<25), l>0DSX0(?5) 
'2, RXO(25) 

COMMON/25/NLCOP2 , ACC2, STEP 

COMMON/26/NLOOP3, ACC3, COTPHI(25), THETA0(50) 

COMMON/?7/DH(50) , HD(25), DDI50J, DRR(50), HH<?5) 

COMMON/28/A, B» C, P 

DIMENSION FOPCEdS), DGDXK251, DGOX2(25)‘ DDDX1(?S), 0DDX7(?5), D 
10(75) 

NAMELI 5T/INPU r/NSHAPL , RADIUS, AZAXiS, BYAXIS, N, I3ETA, DELTAH, AF 

1, YSTARV, ZSIARV, GVORT, D, GS, ICASE, LINC, XO, XF , DX, XCF, XCFs 

2, SPANF, DELTA, DELTAF, EDELTF , ACL, ACCl, ACC2, ACC3 , NLOOPl, ML 
300P2, NLOOP^, LPRES, NPRE, PRPRIT, NP , PDX » IRESTA, LPUNCH, TI^F, 
AMETHOn, force 

20 NW= INLIST (INPUT, LL) 

IF (NW .FQ, 3HEND1 21, 22 

21 PRINT 19 

call PESTAT (LPUMCH, 1) 

CALL CONST 
MP = 0 

IF (LPPES .FQ. 0} GO TO 69 
KC0NT=0 
KDEL=DX/POX 
6P STEP=,001 

ALPHAa=ALPHA/RAU 

XSPF=XCF 

xFPF=ycr 

KST0P=(Xr~Xf')/DX+2. 

IF (LINC .ro, 01 5F, 59 

53 li- (ICASF .rC. 1) 13, 14 
1 3 SPANF = ( XCF-XCFS 1 ^<-DELTAE +XCFS»DFLTA 
GO, TO 7? 

14 SPANF=( XCF-XCFS l^^DFLTAF 
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A, r, P) 


?i t-ALL f IT (X‘>» XCTj ^PAf>;r , DFLIA, PriTAr- 
GO TO 60 

■po ’,iJAiir- = (xrr-:'-rF p)^^For ltt 

•'/'LL TIT {K">» XCF, L^A^'^ , J., FOF(. TF * Ai Hr C ■> Pi 
^'O ir (iprsTA .ecu 0 ) 72 , 7 ? 

72 X-X 1 

>.xn = xp 

XXK=XO 

CAlL NCO*'oH (O) 

CALL SfTUP 

IF (MSh/.pr .EC:). ^' ) 33. 34 
33 L3YAX3 3 = RA0I0S 

IF CBLTA .<^0. 01 -aO 

33 IF (RADIUS ,EQ. 0) 33, 30 ' 

38 PRINT 51 
GO TO 50 

39 PRINT 52 
GO 70 50 

36 PRINT 53 
GO TO 50 

34 IF (BFTA ,F0. 0.) 45, 46 

45 PRINT 54 
GO TO 50 

46 PRINT 55 

50 IF (1CASE .EQ. 1) 15, 16 

15 XY=X 
IPRIT=2 
GO TO 17 

16 XY=X+14,4 
IPRIT=1 

17 PRINT 1, XY,S, A2AXI5, BY AX IS, YO , ZO , ALPHAA 
CALL COF.IGR (1) 

PRINT 2, Gv'ORT, YSTARV, ZSTARV, RSV, THETSV, YV, IV, RV , THETAV, 
L(H( I ) , 1 = 1 , N2 ) 

DO 25 1=1, N2 

2 5 PRINT I, 6S(I), YS(I), ZS(I), RS(I), THETAS ( I 1 , Dill, Y(I1, Z(I 
1) , R{ I 1 , THFTAI I 1 

call lift 10) 
call useful 

DO 26 1=1, M 
L =2^^1-l 

eSPiXO ( I 1 =COSPHI ( I 1 
DWDSXOI n=DWDSGR( I 1 

26 PXO(I)=R(L) 

DO 27 1=1, N 
SUNK I 1=0. 

SUM?( I) =0. 

SUMIOI I )=COSPHI t I 1 

27 SUM2v ( I )=DvJDSGRI I) 

DGVDX1=GV0RT/X'J 
DGVDX2=FORCF ( 3 1 *DELT AF 
0GVDX=DGVDX1 
DYVDX1=Y5TARV/X0 
DYVDX2=F0RrE{ 1 )*PELTAF 
DYSVDX=DYVDX 1 

' DZVDX1=2ST/\RV/X0 

D2VDX2=FORCE (2 ) »DELTAF 

D2SVDX=D2VDX1 

DO 28 1 = ] , N 

DGDXK 1 1=GS( 2*1 1 /XO 

DG0X2 ( I )=rnRCE ( I+o ) *nrLTAF*-*? 

DDDXl ( 1 )=D(2 •f I 1 /XO 
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Df'DX? ( I l-'^On‘CF{ n D-iiOriTAF 
DGDXf [ )^-f;GPXl ( n 
PF Di>r»x ( 1 ) -nooxi ( I ) 
l'K = l 

GO TO 74 

7 3 CALL '?ESrAT (LPUnChU P) 

LX=.k+] 

74 DO < = <K, KSlOP 

GVORTO=GVORT 
YSVO=YSTAKV 
ZSVO = /.STARV 
RV)=RV 
YVO=VV ’ 

ZVO=7V' 

DO 11 7^] » ^!2 
GSO ( 1 ) --6S ( I) 

RO( I) in 

11 THEITAO ( I ) =THFTA( I > 

DO 32 1=1 » N 
DO( I ) =D ( 2*1 ) 

32 COTPHI ( I ) =COSPHI ( I ) /SINPHI ( I ) 

XO = X 

x=x+oy 

XXX=X 

IP (ICASE *€0. 1) 18» 24 
18 XY=X 

GO TO 37 
?4 XY=X+14»4 

37 GV0RT=GV0PT0 + DGVDX«-DX 
Y5TARV=YSV0+DYSVDX*DX 
ZSTARV^2SV0-(-DZSVDX*DX 
DO 30 I=li N 
12 = 2* I 

D{ I2)=D0( ] )+DDDX( I )*DX 
30 GS(I2 ) =GSP ( 12 }+DGDX( I ) is^DX 
CALL RCONSH 11) 

CALL' SETUP 

call loop 3 (ISTOP3) 

IF (ISTOP3.EO. 1) GO to 8 
PRINT 7» XY» S 

PRINT 41,GV0RT» YSTARV> ZSTARV* PSV» 7HETSV» YV» ZV, RVj 
DO 5 I=li N2 

S PRINT 3, 1, OSm» YS(I)» ZSm» RSU), THETAS(I)» D(l), 
1), R{n» THETA! I) 

CALL LIFT (1 ) 

IF (LPRFS cTQ, 0) 42» 63 

63 IF (MP .FO. 0) 64» 43 

64 IF ( XY-PRPR IT ( IPkIT )+DX*.OUl ) 42» 63» 65 
66 DELTAX=OX 

DX=PDX 
MP = 1 

43 K=K-1 

KC0NT = 1CC0NT+1 

IF (KCONT/KDEL .LT. 1} GO TO 66 

K=K + 1 

KC0NT=0 

66 IF (ABS! (XYfDX)-PRPRIT( IPRIT) ) ,LF. .001) 47, 48 

47 call PREPRS (1) 

GO TO 42 

48 IF ( ABS CXY-PKPRI TUPRI T ) ) .LE. .001) 49, 42 

49 CALL PREPRS (2) 

DX=DFLTAX 


THEl AV 
Y! I ) . Z( i 






Mp-n 

li 7 DO ]? 1 = 1, N 

SU'li ( I { I ) 

{ I ) =ssuM? ( I ) 

5,U’'UD( [ ) ■-■ .OMZ 1 { [ ) 

] 2 I } = ' IJI'21 { I ) 

FF {MF'TH'^'' .F(% :i G-' FO 67 
DGVDK= ( G'/-pr~0\/OR7 0 ) /DX 
Di'.SVDX - ( YblARV-^GV' ) /DX 
D/SVDX=-{ZSTARV-X6V0 ) / DX* 

DO 6f> 1=^1 ♦ N 
i 

DGDX t I) = ( G5(L )-(jS 0 ( U ) /DX 
6a DDDXC I ) = ( D(L) -D- { I ] 1 /DX 
GO TO 61 

67 IF FX .LF. XCF) 66. 57 

56 XD= CX~XXO ) /( XCF-XXn ) ' 

DGVDX = DG\/D>^1+XD^ ( OGVDX P-DGVDX 1 ) 

DYSV'DX=DYv/rxH X'J^f ( DYVDX2 -DYV/DX 1 ) 

DZSVDX^DZVDXl eXD> ( DZ VD> 2-D/.VDX 1 ) 

DO 31 I=l» N 

DGDX ( 1 ) =uGDXl ( I '+XD^(- ( DGDX2 ( I )-OGDXl (1)1 
31 DDDXF I )=i)DDXl( 1 HXD*(DDDX2f I )~DDDX1 (111 
GO TO 61 

57 OGVOX^DCV/DX? 

DYS\OX = OYV/DX2 
DZSVDX = DZ\/DX2 
DO 62 I=l» N 
D6DX( r 1 =DGDX2( I ) 

62 DDDX ( I) =0DDX2 ( I ) 

61 CALL RESTAT (LPUNCH, 3) , 

29 CONTINUE 
GO TO 20 

22 IF (NW «E0. 4HST0P) 97, 98 

97 STOP 

98 PRINT 6 

99 NW=INLIST( INPUT, LL ) 

IF (NW .EO. 3HEMD) 20, 99 

8 print 9 
GO TO 20 

1 format (IHO, 23X, 2HX=F10.5» 5X, 2HS=F]0.5, 5X, 7HAZAX 7 6= FI n . 5 , SX 
1, 7HBYAXI S=F10.5//A2X » 3HY0=E12.5» 5X, 3HZO=E12. 5//4AX . 16HANGLF 0 
2F ATTACK = FU).5» IX, 6HDE0REE//) 

2 FORMAT (IHj, 51>, 6HGV0RT=E12.5//20X , 7HYSTARV=L1 2 « 5 , ?X , 7HZSTARV 
1=E12.5, ?X» 4HRoV=E12.5, 2X, 8HTHET AS V=L 12 . 5/ / 2 5X , 3HYV=E12.5» 2XT 

2 3HZV = E12.5, 2X, 3HRV = E12.5, 2X, 7H THETAV = E12 . 5/ / 1 7X , 2HH = 6(E12,5, 

3 ZXl/Zl'^X, 6(E12.5, 2X)//2X, IHI, 5X , 2HGS, 9\ , 2HYS , 9X, 2HZS, 

4, 2HRS, 7X, 6HTHETAS, 8X, IHD, lOX, IHY, lOX, IHZ, lOX, IHR , 7X , 5 
5HTHFTA) 

3 format (IHJ, 12, 10 (IX, EJ 0,3)1 

6 FORMAT (IHI, 20X , ?6U TtilS CASE WILL 8E IGNORED) 

7 FORMAT (1H^,//36X, ■ IHE RESULTS AT X='F10.5, IX, ‘AND S='Fl0.5, IX, 
1, SHARE) 

9 FORMAT (IHO, 73HTHE PROGRAM IS STQPED DUE TO THE NOMCOf'IVFRGF- NC E IN 
] ONE OF THE THREE LOOPS) 

19 FORMAT (1H1, 48X , 24HM0NCON ICAL SMITH PROGRAM) 

41 FORMAT (IHO, 51X, 6HGV0R T = E 1 2 . 5//20X , THY STARV = L"1 2 , 5 » ?> , 7H7STARV 
1==E12.5, ?X, 4HRSV=n2«5, ?X, PHTHI T A S V=E 1 2 » 5/ /2 '>X , 3HYV = El?,5, 2 > , 
2 3H/V-F12.S, 2X, '3HRV = ri2.5, ?<, Ml I HF I AV = F1 2 ♦ 5/ / ?X , IHI, 5x, 2 hc' 
i, 9X. ?HYS, 9X, 2HZ.S, 9> , 2HRS, 7X, 6HTHETAS, 8X, IHD, lOX, IHY. 1 
A.^X, 1HZ» UjX, lH(ii 7X» 5HTHL1A) 

51 FORMAT (3H^, 52X, I6HFLATL PLATE CASE//) 




surj'^oin 1 Nt CONST 


c THIS SUbROUTlNL CAlOILATLS SOME LJStFIJLL C0NSTATNT3. 

COMMON/1 /RAIj'iUS-? AiAXIS, bYAXIS- 5» BETA 

COM.'KJN/^ //>Li'hA\» Ar.o DELTAS OrLTAF» SPmNF» XS^F* TAMAF ♦ IrINAl'. SINA 
1?> COS/-V.?, COSA3* COSHT 

COMMON /h/P i , TVOPI , PIOt>» 1HPI02» RAD, N« NZt N?M1, M?P1 
COMMON/A/GVOPT , GVORTO, bS(SO), GSO(SO), PS(3D, THETAS. SO), YS(SO 
1), 2S(S«1), D(S1), H(52), R(5C1, R0(5G), THETA (50), Yi501, Zf^O)* R 
2SV, THF.TSV, YSTARV, ZSTARV, YV, YVO, 2V, ZYO, RV, RVO , THETAV 
, COMMGn/ZO/XEPF , FDELT) 

’ COMMON /21 /DEL TAH( 50) 

COMWON/27/DH ( , HD (25), DD(.5('), DRR(50), HH(25) 

N2=2^N 

N2Ml=N?-3 

N2P1=N2+1 

PI=3.1 A 15926 53589 793 

TW0PI=2.-^‘PI 

PI02 = P 1/2. 

THpI02=PI+PI0? 

RAD=PI/380, 

EDELTF^TAN ( EDELTF«RAD ) 

D£LTA=TAN ( DELTA^RAD) 

DELTAF = TAN ( DELTAF*RAD ) 

TANAF=AE^-DELTA 
ALPHA=ASIN( TANAF) 

SINA1=SIN (ALPHA) 

5INA?=5IN(ALPHA)«*2 
C0SA2 = C0S ( ALPHA) 5^*2 
C05A3 = C0S (.ALPHA) **3 
C0SBT=C0S(BFTA^RAD] 

H(2 )=DELTAH( 1 ) »RAD 
DO 23 1=2, N 
I2=2«I 

‘'23 H( I2)=H( 12-7 )+DELTAH{ I )*RAD 
H( 1 )=H( 2 ) /2. 

DO 33 1=3, M2M1, 2 

33 H( I )=(H( I-rD+HU-l ) )/?. 

H(N2P1) =H(N2 ) 

H(N2+2)=H(N2) 

DH( 1 ) =H( 2 ) 

DO 34 1=2, N2 

34 DH( I )=H( I+l )“H( I-l ) 

NM1=N-1 

HD( 1)=0. 

HD(2 ) = CH(6 )+HC4) -H( 2 ) ) /4. 

HD(N) = ( 3»*H (N2)-2»*H(N2-2 ) -H(N2-4 ) ) /4, 

DO 35 1=3, NMl 
I?=2«I 

35 HD( I) " (H( 12 + 2 )+H( I 2 ) -H ( 1 2-2 ) -H ( 1 2 -4 ) ) /4. 

HH(1 )=H(A ) 

DO 16 1=?, N 

36 HH( n=H( 2«-I+2)-H(?*I-7) 

RETURN 

END 


.1 Q 



rORMAT 

MMu, 

/i ' 1 y , 

PURPA T 

( 11 iw . 

/(''X , 

r OR.v.A 1 

( IIIP» 

3^X , 

roRi' ,AT 

(111/. 

3DX , 

END 




^'■'tic II'CllLAr< BODY WITH STRAIbHT STRAKfIS//) 

’’•'I iC I RCULA'-’ BODY Vv'ITH DIHr(>RAL MnA*<rS//) 
Wlj, LL iPNCAL BODY WITH SIPAIGDT STPAKrS//) 
OYMriL [PI JCaL B(’DY Vyl m DIHCDRAL STRAKri>//) 



SUHI'Oin iNr PfiiTAT (L« M) 


C 
C 

c 
c, 
c 
c 
c 

CdMI'-'CM/ o/DX» XO, XXX » xxn 

COMMON/5/PU TWOPr» '^I02» THPI02» RAD, R, NR, M2M1, R2P1 
CO.'1MOR/6/GVORr , GVORTO, GS(50), GSO(^jn), ’’5(51) » THETAG(bO). YG ( ''>r 
1), ZS(50), D(51). H(Si2), R(50), R0(50)s. THETA(‘?0)« YC30), R 

?SV, TMFTSV, YSTARV, ZSTARVj -YV , YVO- 2V, ?VO, RV, RVOo THtTAV 
COMKON/7XRO, THETAO, YO , ZO 

COM.’'iO.\!/ll/Dl-'DSGR{2b ) * DWDSGM25), E<25)» COSPHM25), A 

1850(2“^: 

COMMON/12 /X, NPPE, I PR IT 
C0MM0N/17/RADISS, ZaXISS, YAXISS 
COMMON/ 19/DBDX, DPDDX 

COMMON/22/DGVDX , UYbVDX » DZSVDXi DD0X(25), DGDX{25)» K, XY» TIMf, 
IPRPPI T ( 10 ) 

C0MM0N/2A/NL00P1 , ACCl, SUMK25), SUM10(25), SUMil(2^), 6SUMx(25), 

1 SUM2(25)» SUM?0(25)» SUH2K25), SSUM2{2b)» CSPIX0(25). OWDSXO(25J 
2, RX0(25) 

COMMON/29/SUM3, SUMAtJ, SUM31, SSUM3 , SUM4 , SUM40, SUM41 , 35tJ.v4, 3U 
1M5, SUMfjO, SUM?1, SSUM5, SUM6 , SUM60, SUM61, 5SUM6 

NAMFLIST/GOO/K , IPRIT. 6V0R1 , YSTARV, ZSTARV, RV, YV, 2V, G5-, R, 7 
IHETA, D, DGVDX, DYSVDX, DZ3VDX » DGDX » DDDX, X» DX , RADISS, ZAXISS, 

2 YAXISS, XXO» RXU, CSPIXO, DWDSXO» COSPHI, SINPHI, 5UM10, SJMl , SU 
3M20, SU.M2» SUM30, SUiM3» SLJM40, SUM4 , SUM30» SUMS, SUM60, SUM6 , PC, 
4 THETAO, YO, ZO, DBDX , DPDDX 

IF (L .EO. 0) 1, 2 

1 RETURN 

2 GO TO ( 3 , 4 , 5 ) M 

3 CALL STARTIME 
LTIMEO=0 
LTIME=T IME*60000 
RETURN 

4 RTAn 9A, K, IPPIT 
READ 9'=', 

IGVORT, YSTARV, ZSTARV, RV, YV, 2V, (GSd)- 1 = 1, N2], (R ( I ) , 1 = 1, N 
22), (THETA(I), 1=1, N2), (0(1), 1 = 1, N2 ) , DGVDX, DYSVDX, DZSVDX , U 
3DGDX(I), I~l, N), (DDDX(I), 1 = 1, N)» X, DX, RADISS, ZAXiS-S, YAXISS 
4, XXO, (RXO(I), 1 = 1, N), (CSPIXOd), 1=1, N), (DwnsXDU), 1 = 

5] , N), tCOSPHl(I), 1=1, N), (SINDHim, 1=1, M), {SUMIO(I), 1=1, N 

6) , (SUMKI), 1=1, N), (SUM20m, 1=1, N), ( SUM2 U ) , 1=1, N), SUM’O 
7, SUM3, SUf‘i40, SUM4, SUM50, SUMS, SUM60, SUMS, PO . THETAO, YU, ZO, 
0 DBDX, DRDDX 

CALL OUTLIST (SOO) 

RETURN 

5 LTIME1=LAPSTIME ( DUMMY ) 

IF ( (XY+2.»DX>'PRPRn (IPPIT) .LE. .01) 44, 45 

44 ir (LTIME -4.*LTIME1 +3.^'LTlMrO“300nn )7n, 71, 7] 

45 IF (LTIMF-?"«LTIMFl + l TTMFO-30OnO) 70, 71, 71 
71 LTtMf o=lt run 

RETURN 

70 PUNCH 06, K. IPRIT 
PUNCH 95, 


NOT REPRODUCIBLE 


THIS SUP'*OUTjNE pi ads A/JD PUNCHES NECFRSSARY DATA FOR RTSImPTINC. 
THF PROGPn'>ri NCNKOnSM. 

1=0, 1 UE_PR('GI’AM NOMCc'iNSM DOFS NOT NFED TO RESTART. 

L--1, THL PRwGPAN NONCONSM L'ILI BL PLSIARTLD AT aNY TIME. 

M--1 * PRLPaRLs t-UK KE'S 1 APTING. 

M-2, RlADo AN(; CliECkLS IHE REST/\RTING DATA. 

M=3, PUMCHfcS OUT TUL NECERSSAPY DATA FOR PFsrARTiNG. 
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lOVORl * Y.ST<\RV, -ili,7ARV, RV > YV» {GS(D. 1 = ], iN2), (R(I), 

?2), nuriAU)? 1^.1, N?).» (DMl* N2 1» DC-vnx> OY'.VlJX, D/i,VI. 

'^D{.r>x(lM 1 = 1^ '7), (i)i)DX(I), !=■]. N), X» DX» RADlOSs Z-AXISA, Y/. 
iyy XXO. IRX 1(1)9 I--1. Ml» (Ci>PtXO(I)> 1=1. N), ( Dv'i)' I ) . 
'll. N). (rObPHlil), r = l-. f't)9 (SIMP^HI(I)» I=l9 N), (AUMin(j), 

61 » ISU'Utn, i-1, N)» (GUM?J(I)9 1=1 9 N), (SUM7(1), 1=1. 

7, SU’-H, AL'MAr,, SUMA , r.UM5» SUN'60, AU-16, RO. THFlAO, YO i 

8 DBDX» DRDDX 

9A FORMAT (1X» 3(L23.15-, 2X1) ^ 

FORMAT (3b23.]‘n 
96 FORMAT (21 1^) 

STOP 

END 


I 9 M 

V , |J 

.XIO'- 

7 9 . 

1 9 .0 

ZTS 
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SUBROUTLNf FIT (XO. XCF, SF > DELTA, DELTAF, A» B» C, P) 

THIS Sin’PIHJTIHF Firs TIIL CURVE OF MONCONICAL SECTION BY A 
nURO OROFP POLYNOMIAL, 

Sn=xO^Dr LTA 
XO?.=XO»^5"? 

XF2=XCFrf*2 

XF3=XCF**3 

A3=XF3-X03 

B3=XF?-XO? 

'C3=XCF-X0 
F3=SF-Sn 
Al = 3. »XO? XC3-A3 
B1=?,«X0«C3-B3 
Fl^DELI A-WC3-F3 
A2=--3,*B3 
B2=2.*C3 
F2=DFLTAF-DELTA 
DV = A]*B2-A2^fBl 
A=(B2*Fl-BI^fF2)/DV 
B=(A1*F2~A2*F11/DV 
C=<F3-A3 kA-B3!^Q) /C3 
P=S0-X03*A-XU2^8-X0«C 
RETURN 
END 



SURWUTINt: NC0M6H (L) 

C THIS SUB'"’0U1 INi: CALCULAT! S SOME INPUT DATA FOR FACH SECTION (\f K. 

C l. = 0, roR X=XU. L = l, rOR ANY X. 

COMMON/ I /RAO I US » A/aXIS* BYAvXISj S» BETA 
COMMON/ 2 /NSHaPF. 

COMMON/3/DX, XO, X, XXO 

COMMON/A/ALPMA* AL. DtLTA» DLLTAFi SPANF, XSPF', 1AMAF« SINAI. S I NA 
12 » COSA2» COSA3» COSBT 

COMMON/S/PI* TWOPI, PIO?» THPIO?» PAD.. N, N2 » N2M1 , M2P1 . 

COMMON /A /OVOR r , GVORTQ» GSISOI, GS'MSO), RS(SI), THTTASISO], ys(‘>n 
])5 ZS(SO)? D(S15» H(S2), R(SO), R0(5O), THETACSOl, 2 ( 'O ) , r> 

25V, THFTSV, YSfARV, ZSTARV, YVO, 2V, ZVO, RV, RVO, THlTAv- 

COMMON/9/EPS 
COMMON/ 10/L I NC 

COMMON/17/RADISS, ZAXISS, YAXISS 
COMMON /19/DBDX, DRDDX 
COMMON /20/XEPF, EDELTF 
COMMON /2 3/ I CASE 
COMMON /2R /A, B, C, P 

IF (L .EQ. OJ 1, 2 
1 Xl=l. /DELTA 
RAT 10=XX0/X1 
EPS=0, 

IF INSHAPE ,rO. 0) 4, 5 
4 RADIUS=RADIUS«RAT 10 
% ORDDX=RADIUS/XXO 
RADISS = RADIUS 
GO TO 6 

l5, AZAXIS-AZAXI S*RATIO 
_;?• ZAxrSS=AZAXI s 

BYAXlS = BYAXIS-*fPATIO 
DBOX=BYAX IS/XXO 
V; YAXISS = eYAXlS 
^6 YSTARV=YSTAR /*RATIO 
*/- 2STARV=ZSTARV-K-RATI0 
?* DO 3 1=2, N2 , 2 

65 { I >=G5I 1 J^RATIOif-IDELTA) 

3 DU )=D[ I )*RATIO 

GVORT-=GVORT*RATJO^<-DLLTA 

S=RA1 10 

RETURN 

? IF (LINC »F0, 0) 19, ?i 

19 IF (XeLh. XSPF) 10, 11 

10 S=A«X**3+B’-^X^i'*2 + C*X+P 
9 RATIO=X/XXv 

IF CNSHAPF .LQ, 0) 7, 8 
7 IF (ICASE ,FO. 11 12, 13 

12 RADI US=RADISS*RATIO 
GO TO 20 

13 RADIUS=S*?,/3, 

DRDDX^^Z.wiB.-^^A-s-X »*2+2 . *B»'-X+C ) / 3 » 

20 RETURN 

11 S=SPANF+(X-XSPF) ^<-DELTAF 
RATIO=X/XXU 

IF (NSHAPr .1,0. Oj 14, 8 

14 IF (ICASE .FO, 1) IS, 16 

15 RAOIUS = PAnrsS^t-RATTO 
GO TO 20 

16 RADTUS=-S*2,/3. a 
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Gn TO ?f> 

A2AXJ'.=/-AX1SS^RA1 10 

BY AX i b = YAX 1 GSx^KAT lO 
P^ 1 ’JRM 

PAT TO^-X/XXv' 

RAO T U.S = f’ AD I S SifR A T I 0 
S-v-wDPLTA 

TF (X .LF. XFPF) 17, 18 
FPS = 3e * A-« X P+2 « ^B+X ^ C 

PFTURN 
r.P8=^FDELTf- 
^ RETURN 
-FND 



SlIBROUTIMC 5f:TUP 

C THIS SUt3ROUTlNL- 6LTS UP ALL C OWSTAI'JT S FOR THF T i' ANSFORM a T I ON 

COilMON/ 3 /KaH I'Jb » AZAXIb^ BYAXIP^ St nriA 

CO'1MOK/P/i\SHmO[ 

COMPON/5/P I ♦ 7-(OPl, PIO?, 7HPIO?, RAD, ^!, M2- 
Cn*'1MOM/7/PO , THETAO, YO, 70 
CO^^t'IO.M/a/SPRSQS, SMR5C3S 

COMPON/13/ LA, CB» CBMA, C)jMA,S£>» L FI A2 , CPFO 
C0MM0N/14/AZ , BY, A2S0, bY^O, XIZ, C A2 , RMA, bSOAS, RSO 
COViWON/15/ A, bs- ROOTA, SI.'IANO, COjANG, BtTAH 
COf/.MON/30/TLlFT, UPPER, XLOlv'Fi^, 

COMPLEX CA, CB, CBMA , CE-'fYASQ, CFETAZ- CR50 

IF (NSHAPF .EQ. 1) 1, 2 
1 AZ=AZAXIS 
CA=AZ 
BY = E3YAX IS 
CB = BY 

A750=AZ*->'2 ‘ 

BYOO^PY^i-KZ 
8ETAR=BETA*RA0 
T = 2.-«BETAR 
BMA=BY-AZ 
CBMA=BHA 
.'•V BMASQ = BMA*«? 

't'CBMASO = RMASO 

BSMAS = BY^*?-AZ*-i<-2 
A-C0S(T)-BSMAS/Si«-»2 
\% B=SIN(T ) 

POOTA = SORT(SQRT( A**2+B*-^2 1 ) 

'^.-\hET=ATAN2(o> A) 

■“£MF.(THET .LE. 0.) 3, A 

3 -/thet=thft + twopi 

4 ';;^ANG = .3*THFT 
^SINAMG = SIiM(ARG) 

COSANG=roS( ANG) 

FAcf=PY*R00T4 

X IZ=Sf ( FACT-R-COSARG-AZ^COG IBFTAR ) ) /BRA 
X IZSQ=X IZ»^^2, 

ETAZ=3R < FACT^'-SINANG-AZ*SIN { BF TAR ) ) /BMA 
C 1FTA7 -=FTAZR ( 0, , 1.) 

TLIFT = X IZ R^?+BY5Q 

U P P t R = 3 w R T U E G V A G - H ' 1 A £ T A Z ) R » 7 / ■ ' 1 A S C J + X 1 Z S 0 ) 

XLOV^FP^ -nU'-T { ( P5,^ AG-*I’!’A-vF 7 AZ ) iS'VPrAGQ+X 17 30) 

OO TO 7 

7 RSQ = RADruG*!-7 
CR50=RGQ 

IF {RSO .rO, 0.) 5, 6 
3 RLTA=J. 

6 b[-7AP=BrTA«-RA0 
5‘^RGOG-^^ + RGO/S 
FTAZ = ( G^PS0/G)*SI^1(BETAR) 

C i r TA7-=^f T A7‘f ( O. , 3 , ) 

S^’RGOG-S-RGO/G 
X 1/ = ( G-RGO/G) «C0G(B1 TAR) 

X IZGO = X 1Z*<»2 
TL IFT = XlZ»*7<PGO 

UPPFR='i JRT( ( ?.‘‘RAUUlG-f T A7 ) «*?+XlZSO) 

XLOWf R --GORT ( ( 7. «PADniG+f lAZ )^^«? + X I /GQ) 
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CAu y/rovu (M,, ^ 

p+2')i> V p ) 

"T MF T A0=ATA,\| f / 0/ Y Q ) 

RMtJRM 

TRD 
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;;uRf?(nniwr <y, yS» c5) 

THrs SifnROi/n.'iE Dors ruf transformation rstwfe,n z and zstar 

rci‘'MOM/2/N5'(APF 

c ;miu‘;i/e»/p I , T'A;pj , ptoz, mpjoz* pad, n. nz, nzmi, n?pi 

A /y\z . BiT, A/jO, i’.YSO, XlZ, tiTAZ, OKA, HSMAS, RSO 

IF (NSliAPF .FO. 1) 1, 2 

2 DrNO'M= Y««2-+Z-"-«2 
pe.L=Y-RSO*Y/OrNOf^ 

A I M= Z + R SO’f-Z / DFNO.v-F T A 2 
T0P~2.5<-REL^^A IM 
hot TOH=REL^f'^' 2 -A I .v^cir2-X I Z^-*2 
RS=SORT ( SORT < TOP^ * 2 + BOTTOf/-s^-K-? } ) 
rF (TOP) 8, 12, 8 

12 call AMGDFT CRFL, aim, phi, bottom, COSTRv, SINT'?''! 

GO TO 10 
I T0P = 2.^Y-'^Z 

BOTT0M^Yv< 2-2^-^'-2“UY50+AZS0 
R = SORT f SORT ( TOP*»2+dOTTOM^<-»2 ) ) 

IF (TOP) 4, 5, 4 

5 call ANGDET (Y, Z, phi, BOTTOM, COSTRM, SFNTRH) 

GO TO 6 

4 THET==ATAN2 ( TOP , BOTTOM) 

IF (T!in *LT. 0. )THFT==TH>^T + TWOPI 

PHI=.5-»THET 

COSTRM = COS(PHI ) 

SINTRM=SIN(PHI ) 

IF (COSTRM -DC. 0.) 6, 7 

7 C0STRH=~C0STRM 
STNTRM=-SINTRm 

6 REL=flY-^R*C0STRf1-AZ*Y 

aim=by»r-)^sintrm-az-»z+bma^<-etaz 

T0P=12.*REL«AIM) /BMA-'^-if-Z 

bottoms: (REL^'^Z-A I M**2 ) /BMA**2-X 2 

RS-SORT (SORT ( TOP^^i-2+eOTTOM'f^f-2 1 ) 

TF (TOP) 8, 9, 8 

P call ANGDET (RFL/BMA, AIM/BMA, PHI, BOTTOM, COSTRM, StNTPM) 
GO TO 10 

8 THET = ATAN2 (TOP , BOTTOM) 

IF (THE! ,LT. 0.) THE T=THET+TWOP I 
PHI=. 5«THET 
C0STRM=COS ( PHI ) 

SINTRM^SIN ( PHI ) 

IF (COSTRM .GF, 0. ) 10, 11 
n COSTRM=-COSTRM 
SIMTRM=-SINTPM 
10 YS=:PS«COSTPM 
ZS=RS«SrNTRM 
RE TORN 
END 
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U'’,RO!.n If'T (YS, ZG » 


Y» Z) 


c 


2 

] 

li 


6 

3 

8 


,9 


li 


15 
17 
14 

16 
7 

19 

IB 
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THIS M'BRO'n INC- Dt)f-3'THr T('ANSFO«^ATIO/^l F3F.TW^:£^^ Z'^TAP Af'n Z 

i 

CO'lMnn/Z/NSrfAPF 

LOKM0(>:/5/P I » IVuPI, fMU2» IMP I 02* RAD. I'J » iN 1 , NZl'l. .‘•JZP L 
C OMPiOM / 14 / AZ » R Y > A Z30 , 07 . X ( ? » L T AZ » (5v.A , BSfV> S » R30 

T0'^ = ?.*Y6*ZS 

BOT TOK= YS*-* ?-Z5*«2 + X I Z^«2 

KHO -SORT [SORT ( TOP** 2+ BOT ! OM**? 1 ) 

IF (TOP) 1. ?» 1 

CALL AMOOFT (Y5, ?S» PHI. BOTTOJil. COSTRM. SIHTKM) 

GO TO 3 

THrT=ATANZ (TOi-->, BOTTOM) 

IF (TriFl .LT. 0.) 4. 5 
THF T=THLT+TwnP I 
PHt=.5ilHFT 
COSTRM=COS(PHI ) 

SINTRP. = SIN(PHI ) 

IF (COSTRi'^ .OF. 0. ) ' 3. 6 
C0STRM = -C0STR'1 
SINTRM=“5INTRP 
IF (NSHAPE oEQ. 0) 7, 8 
AL = BMA*Rh?0-'fCOSTRM 
BE=BMA*{RH0*5INTRM+ETAZ ) 

TOP=?.*AL*BE 

B0TT0M=BSMAS**2+AL**2-BE**2 
RHO=SQRT ( SORT ( T0P**?+B0TT0M**2 ) ) 

IF (TOP) 9?1J» 9 

CALL ANGDET ( AL t BE » PHI. BOTTOM. COSTRM. SIMTRM) 

GO TO 1 1 

THET=ATA|vJ2<TnP. BOTTOM) 

IF CTHFT .LT. 0.) 12. 13 
THrT = THFT-»-TWOPI 
PHI=. 5*THET 
COSTRM=COS(PHI ) 

SINTRM = SIR(PHI ) 

YT.FRM = RHO«COSTRM 

ZTERM=RHO*SINTRM 

Y=( AZ*AL+DY*YTERM)/BSMAS 

Z= ( AZ*BE+f3Y*ZTERM ) /BSMAS 

IF (Y) 14» 15, 15 

IF (AZ*AL'-BY*YTERM) 16, 17, 17 

IF ( ( (Y/BY1**2+(Z/AZ)**2+.U001 ) .GE. 1.) 16, 14 

Y= ( AZ^-AL-BY^YrERMI/BSMAS 

Z= ( AZ-KRE-BY*ZTERM)/BSMAS 

CONTINUE 

RETURN 

A=phO*C'''STRM 

L3 = RHO*SFNTR‘1*ETA7 

IF (RGO) 18, 19, 18 

y~-A 

Z-B 

GO TO 20 
T0P=2.*A»n 

BOT TOM- A*!t2-n«*?+4.*RS0 

PH0 = 50RT(S0RT( T0P*■>'•2^ ROT TOM**? ) ) 

IF (TOP) 21, 22, 23 

CALL ANGDET (A. 13, PHI, FjCfTOM, COSTRM, SlNThM) 

GO TO 23 

THET=ATAN2 (TOP , BOTTO''!) 
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Jf HHL r .|,T. 0 . ) 2 A, 2 ') 
Tin r rif'T-t I WOP I 
?h PHI=.‘'5fTnn 

CO*;! Rf/vC OS ' ) 

SI/m1Rm=SIN (PHi )■ 

2"} YTCPM=RHo»co'^rr?f'> 

/Ti:R'1 = RHO^f^M0TR ^ 

Y = , 5« ( A^ Ylf fM ) 

7=,S«^^^^ZTf P^] 

Jl- <A*-YTf:RM) 26, ^7, p/ 

?1 ir (A-YTFRM) ?0, 28, 2fi 
?a IF { ( Y^^^^7^2ff•t< 7 + .U( ()] ) .OF. 
Y= , S * { A-YTLRR J 

Z=,5*(0-ZTCRM) 

20 CONTINUE 
RETURN 
END 


RSO) 20, 26 
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•■A/BROUTINr DFRIVF lYS, Zb, Y, 2. DR? Dll 
C IHIS Si'BRnuT-INl: CALCULATES DZ/D2S 


CO''*MO(']/?/I.'SHAL-'E- 

rOMMON/b/PI, T>OPf. DI02, THPI02« RAP, N, W2 , N2Mi» N2P1 
COvMnNrn/ CA. CB, CBMA, CBMASO* CJFTAZf CPSO 
CO^iMOM/lA/AZ » By, AZSO, SYSQ, XIZ, ETAZ, PiAA, bsmas, rso 
COM f.^nN/16/ROOT 

CO<>^PLEX ROOT, CA, CB, C8MA, COMASQ, CIETAZ, CAPZS, CAPZ, PZDZA, CP 
ISO, FACTOR 

CAP7S = '''''^LX f YS, ZS) 

CAP/=CMPLXfY, Z1 

rr (NShAPE ,t'Q. 0) 1 , ? 

? TOP = Z.^^Y-^Z 

BOTTOM=Y>^2-Z**2-BSMAS 
RHO=SQRTlSQRT( TOP^’^Z+BOTTOM^WZ ) ) 

IF (TOP) 3. 4, 3 

4 CALL ANGDET (Y, Z* PHI, BOTTOM, COSTRi^!, SINTRM) 

GO TO 5 

3 THET=ATAN2 (TOP, BOTTOM) 

IF (THFT ,LT. 0.) 6, 7 

6 THFT=THFT+TWOP I 

7 PHI=.5*THET 
COSTRM=COS(PHI ) 

SIAJTRM=SIN(PHI ) 

IF (COSTRM .GE. 0> 5, 8 

COSTRM=-COSTRM 

SINTRM=-SINTRM 

ROOT=CMPLX (RHO^f-COSTRM, PHO*SIRTRw ) 

DZDZS= ( CBMAS0*CAPZ5*R00T ) / ( ( CB-^ROOT-CA^-CAPZ-CBMAS-Cl ETAZ ) ( CB-)"CAPZ- 
1CA*R00T ) ) 

GO TO 9 

FACTOR=CRSO/CAPZ 

DZDZS^CAPZS/ ( (CAPZ-FACTOR-CI ETA2)*( fU,0, ) +FACTOR7CAP2 ) ) 

DR^REAL (D2DZS) 

DI=A1MAG(DZDZS) 

RETURN 

END 
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fjIlBROUTIMF ANODE.T (Y, Z, PMI ♦ R- COSTRM> SI'^TRM) 


C THIA SlinPOUTINil DETI'F^'MINf-S AliALU OF i f Y+ I Z ) ■*< « RFAL) 5^-i- . AV.WHF N' Ylop 

C 7 KODAlO ZFR*^ 

COMF.On/'j/PI » rWOPI, PI 02 * THPI0?» RAD» N, N2, NZM1» NZPl 

IF ( R 1 ] , 2 , 2 

1 IF (Y) 5* 

3 IF tZ 1 A, 2, 3 

3 PHI =P I 02 
C05T'?M=0. 

SlMTRM^l . 

PE TURN 

? PHI=0. 

Cn.ATRM-l . 

5lNTRM=n, 

RFTHRM 

4 PHT=THPln2 
C0STRM=O. 

SIN1TRM=^-1 • 

RETURN 

FND 
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iUhlROUIINf COriGR (M) 


^[^lb ';uRRr)UT [RF calculatls all coRr i gurai ions Oh sums and 

C Tlir ISOLATLD VORTfcX 


C A = i , 

r 


GJvr YSV, 
CAl CUTh 


/sv, n(Piv)» o:piv), sv, osipjv) 
ALL Tllh: OTHI RS 


C M = 0* 


calcolatf g(piv) only 


CO.VWON/3/OX, XO» X.. XXO 

COMMON/S/PI, rwOPI, p]02» THPIO?^ RAD i N» NZ? N2M1* N2PJ 
CfjN,MON/6/OVnRT» GVQRTO» GS(ST)? GSO(SO(» RS(Sl)* THC TAS ( "' 0 ) * YSISO 
l)j ZS(SO), R(50) ,-RO(SO) , THETAtSO), Z ( *^0 ) i !< 

2SV, THFTSV* YSTARV, Z5TARV» YV» YVSt ZV, ZVOi RV, RVO» TriFTA\' 
COY'^OM/7/RO, THETAP, YO, ZO 

C0MM0N/27/DH < SO ) * HD(2?)» DD(SO)» ORR(SO)» HH{?5) 

NM1=N-1 

IF ( ,EQ. 1 n » 2 
1 RSV = SQR T ( YSTARV*»2+2St ARV«’!-2 ) 

THET5V=ATAN (Z5TARV/YSTARV) 

call YZCOMP (YSTARVj ZSTARV, YV4 ZV) 

RV=SQRT ( YV^^^2+ZV*-K2 ) 

THFTAV=ATAN(ZV/YVI 
D( 1 ) = ( RSV+D( 2 ) ) /2 . 

ARG=THFTSV+H( L ) 

YS( 1 ) =YSTARV~D ( 1 )*COS ( ARO) 

ZSm=ZSTARV-D(l )^SIN(ARG) 
i RS( 1)=^S0RT(YS( 1 )«-s-2+ZS( 1 )^^2) 

THETAS I 1 ) =ATAN (ZS ( 1 ) /YS( 1 ) ) 

call YZCONt^ ,<YStl)» ZS ( T ) , YIl), Z{D) 

R( 1 1=S0RT(Y( 1 )**2 + Z( 1 )*-»2) 

' THETA ( 1 )=ATAN{Z { 1 ) /Y ( 1 ) ) 

1 . DO 12 I-S, N2R1 » 2 

12 D(I)= (D( I+i )+0f I-l) )/2. 

DO 13 1=29 N? 

APG = THFTSV+H ( I ) 

Y5(‘l )=YSTARV-D( I )-»-COS(ARG) 

ZS( 1 )=ZSTARV-D( I )*5IN(ARG> 

RS( I )=SQRT (YSI I )ft*2+ZS( I )^^*2) 

THETAS ( I )=ATAN(ZS( I )/YS( I ) ) 

CALL YZCOMP { YS ( I ) » ZS ( I) » Yd), Z(D) 

R( I )=SORT< Y( I )**2 + Z( I )**2) 

13 THETM I )=ATAN(Z( I )/Y( I) ) 

23 D(N2P1)=D{N?) 

R (N2PI ) =P (W2 ) 

DD( 1 )=D{2)-R5V 
DRR ( 1 ) =R( 2) -RO 
DO 15 I =2 9 N2 
IPl=l+l 
IM1=I-1 

DD( I) =D( I PI ) -D( IMl ) 

15 DPR ( I ) =R ( IPl )-R( IMl ) 

GS( N2-1 )= ( DH (N2-1 )+DH(N2-’ ) ) ^ ( GS ( N? ) +GS ( N 2- 2 ) )/( 2 . * ( DH I N2 -1 ) +PH ( N? 
1-3) ) ) 

DO 19 1=29 NNl 
I2=2«(M-I+1) 

19 GS( ! 2-) ) = ( ( DH( 12-1 ) +DH [ 12 + 1 ) )-»GS{ 12 )-DH( I 2-1 ) *GS ( I 2 + 3 ) ) Zf)H( I 2-f 1 ) 
GSn ) = ( H< A } •'GSI? )-2.*H ( 2 ) ^->^?*-GS ( 3 )/(H(2)+H(4)))/( ? .>DH ( ^') ) 

return 
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? (,S( ? ) s f ? . XHH r ^ { 1 ) ^ p ( 3 ) / tHf 7 )+H ( ^( ) ) ) /!! ( A I 

[•>0 in r-7» ^ 

: ?=? » ’ 

iS U->( ]/> = (U'H I 7-1 1 ■it'-.nc IPM l+LKU J2+n^GG( I 7-1 ) 1 / (DH( ) +[>m 1 ^ i ) ) 

C.S( N,i) - ( (p. + piKN?-! H-OIUN7-3 ) ! ’'^GS(N2-I l-DH( \7-l ; »^(1S [ (i?-3 ) 1 7 ( Ot I [ ^ 7 - 
1 1 )4DH ( M7-3 > I 
rUFTUR'J 
END 
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surkgi,'tim: sojrcc (v, dii sn* sd 

C THIG Sl'bROUTlMf CAl. CUL /^T h R THT SOURCt ITRM IM COiXPLL X VPtOCIlr 

COMMON /■’/flOHAPf. 

CnMMO^J/l^/ CA, CH» CRf'A» CHMAGOt CIFTAZ^ CPV) 

COMMON /I d/root 
COMMON/ 19 'DaOX , DPDOX 

COMPLLX ROOT, DZ02S, S« CA , CB", CilMA, CBMASQ, C1LTA2, CRSO, CAPZ 
ICDRDDX, CDBDX 

CAP2=C'>'PI X (Y, Z) 

D702S=rMrJLX(DR» DM 
IF {^'SHAPT .fc-Q. 0) 1, 2 
?. CDBDX=DBDX’*{ 1. ,0o ) 

S=CA*CDBDX^^D2DZS/R00T 
GO TO 

1 CDRDDX = DRDDX*{ 1, ,0. ) 

S=C50PT ( CRSO ) »CDRDOX*D2DZ5/CAPZ 
3 SR=REAl(S) 

S r = AIMAC7(5) 

R FT URN 
END 
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SUOROUTlMf-- PPSLH (VS, 2S, Y, DR, DI, Fp . £j , 


THIS iiUf 


BKOUTlMf CALCULATLS liPS = “. ( i D22 /DZS2 ) / { Dr /OZ S ) 

COMmdM/ z/r- SHAPE 

C.OYK’.ON/ 13/ CA» CB, COt'IA, CBMASQ, CIETAZ, CRSD 
CO'-1i'‘r'M/J6/PPOT ■ 


complex CA, f.B, CBPA, CRVASO, WFTAZ, 
IZS, EPS. CRSO» Tl, T2» DI , 02 


RCn T , 


V, .*/, CA»ZS, CAP/. 


rzn 


2 


1 


3 


CAPZS=:CMPLX{ Y6, ZS) 

CAPZ=CMPLX [ Y , 2} 
D2DZS=:CMPLX(DR, DI) 

IF (HSHAPE .to. 0) I, 2 

V=CB»R00T-CA^-CAPZ-CIETAZ*CRMA 

W=CB*CAPZ-CA*RnOT 


EPS=(~. S,C;. ) « ( (1 
I /vn /ROOT) 

GO TO 3 
T 1=CR5Q/CAP2 
T2=TJ /capz 


. .0. ) /CAPZS+DZDZS*(CAP2/P00T-M/V~fCBxf?00T'-CA^.CAPZ) 


Dl=:CAPZ-T]-r]ErAZ 
D2=(l,,0, l+T? 

FPS=( .3,0. ( (D2/Dl-(?, ,0. 

fp^pFAL f FPS ) 

EI = AIMAG(EPS ) 


» (T2/CAPZ )/D2 )sDZDZS-n ,,n. ) /CAPZS) 


RETURN 

END 
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.SliRpniiTIM- DW'IMC 7'',. plD7^P. 

C THIS ^iiPROHnivir f/'LCULATI- S THr TERP 0'/'(IRC)/n?'- 

COf'i/i'/,M/l/KAniUS. A7AXLS. iiYAXIS, S- L3f TA 

COMrOP/'>/RI » P102* THPIC'^, RAD, N, N2 , , .N2P1 

c OMvoft/8 / ,:.PRbC'S , dmr: JS 

CCViVOM/R/T°S 

CO'IPLpy CZS, CA, CR, CD. CE , CF, CT, CDJDZS, CFPS, CTT, CPI?- C^I 

TT=ACOS( 2. -^R ADJUSTS / ( S^--*2+R AD I US^< ) ) 

IF ((YS .EH. C.) .AND. (ZS .FO. 0.)) 6, 7 

6 D1DZSR=0. 

D IDZSI=“EPS* ( TT/P I ] 

RETURN 

7 CZS=^YPLX{YS, ZS) 

CD=CMPLX(0., SNRSOS) 

CERS=CMPLX (FPS, 0.) 

CPI?=CMPLX(TWOPJ , 0.) 

CTT^CRPLX , TT) 

CPI=Cf''PLX(PI , n.) 

CE=(?.«RADlUS*a., 0.])*rzs 
TOP = 2.i^ZS«-YS 

R'^TTOM = YS^<->*^?~ZS*’'‘2 + SPPSOS»’"2 
RS = 50 RT (SORT ( TOPk* 2 +BOT 2 ) } 

IF (TOP) 1» 2^ 1 

2v CALL ANGDET (YS, 7S, PHI, BOTTOM, COSTRM. Slf'TRf/) 

TO 3 

r|;<THET=ATAN2(TOP, BOTTOM) 

IF (THFT .LF. 0.) THET =^THET + TWOP 1 
PHI = .5«-THFT 
COSTRM=COS(PHI } 

.V,S INTRW=S1M(phI ) 

^F-IF (COSTPM .GE. 0.) 3, 5 
5 ‘^JCOSTRM^-COSTPM 
U'SINTRM = -5INTRM 

3 YCT = CMPLX( RS^^COSTRM, RS»SINTRM) 

Cr = CD*CT 
CA=CLOG(CF+rF ) 
rB=rLoo(CP-rF) 

AAR=RFAL <CA ) 

AA I =^AIRAG(CA ) 

IF (AAI ,LF. O. ) AAI=AAI+TWOPI 
CA=CPPLX( AAR , AAI) 

RBR = RF AL( Cfi ) 

BBI =^AIMAG(CB ) 

IF (HHl .LE. 0.) BRI=^BBI+TW0P1 
CB-^^CMPi X (PRf5, BFI) 

C0IDZS= CFPS^[CZS*(CA-CB) / ( CPI 2»CT )-CTT/CPI- (0. , ].)*((]., 

1/CT )/(?., ■).)) 

DIDZSR-RFAL(CDIDZS) 

0 1 DZS I = A 1 mag ( CD I DZS ) 

RETURN 

FND 
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r>liRPOUF INF USFFUL 


IMI.S SUl'POUriNf CALCULATtS SOMt USFFUL EXPRESSIONS FOR 
TilE (OOPl AMD THE LOOPS 

COM^;ON/S/DX . XO, X» xxo 

COM;'’ON//i /alpha , AF , OFlIA, DLLTAF, SPANF, XSPF , TANAF, SINAI, SIN'- 
i;s COSA?, COSA3-. COSRT 

COM-'''ON/S/PI , TWOPI, t>102, THPIOP.f J^AD, N» N2> N2v.i, n.^pi 
COMMON/S/GVORT» GVORIO, GS(50), GS0(50), RS(51). THETAS(S01-> 

1)» ZSCH'}, D(bl), H(52). R(bO), R0(50), TMETAISO), Y(50), 

2SV, THFTSVj YSTARV, ZSTARV,- YV. YV0» ZV, ZVO, RV, RVO, iHETAv- 
COMRON/7/RO, THE T AO. YO , ZO 
COMMON/ 10/L INC 

COMMnN/ll/DUD5GR(2S) , DV.'DSG I (>S ) , F{?'^), COSPHI(?F). A 

l'SSO(2b', 

C0-MM0N/27/0H ( SO) ♦ HD(?S), DD(SO)» DRRCSO). HH(2b) 

DIMENSION UR(25)» DI(2S), DWDZSR(25)» DWDZSI(25) 

DO 3 1 = 1 » N 
- L=2«l-1 

CALL DERIVE ( YS ( L ) , ZS(L),Y(L), ZfL)» DR ( I ) , DI(in 

3 ABSD( I )=] . /SORT! DR { I )*^^2+Dl (I )**2) 

DO 1 1=1, N • 

I2M1 = 7^"I-1 

CALL SOURCE {Y(I2MD, Z(I2^M1), DR (I ) » DI(I), SR, SI) 

-,A-f (LINC .E.O« 0) 6, 7- 

6 pIDZSR=0. 

.■i|’I-DZ5l = 0, 

■G,p‘ TO 4 

7 ©All DWINC (YSTIZMI), ZS(I2MD, DIDZSR, DIOZSI) 

4 ;Y*SMYSV = YS< I 2MU -YSTARV 
.Y^SPYSV=Y5( I2MD+YSTARV 
•gSMZSV = ZS( I2M1')-ZSTARV 

'i|)ENV= ( Y5MYSV**2+Z5MZSV**2 ) * ( YSP YSV**2+ZSMZSV* *2 ) 

'SHijMR=0. 

.IDmi=o, 

■D'O 2 K=1 , N 
K2=2*K 

YSMYS=Y5( I2M1)-YS(K2) 

YSPYS=YS( I2MI )+YS(K2) 

ZSMZS=ZS{ I2M1 )“ZS{K2) 

YDENS= { YSMYS**2+ZSMZS**2 ) * ( YSPYS**2+ZSMZS**2 ) 

SUMR=SUMR+GS ( K2 ) *HH ( K ) * YS { K? ) *ZSMZS^ I YSMYS+YSPYS) / ( TWOPI*YDFNS ) 

2 SUMI = SUMI+GS (K2 > *hH ( K ) *YS ( K? ) * ( YSMY S«^YSP YS--ZSMZS**2 ) / ( TH'OP I«YDEMS> 
DWDZSRt I )=-GVORT *YSTARv*ZSMZSV* { YSMYSV + YSP YSV > / ( P I -s^-YDENV ) -SUMR-+- 

ISP+DIDZSR 

1 DWDZSII I)=-(GVORT «-YSTARV>^ C YSMYSV* YSPYSV-ZSMZSV**2 ) / ( P I *YDI NV ) +S 

lUMI +TANAF-SI -DIDZSI ) 

DO 5 1=1, N 
L=2»l-1 

L ( I ) .-SQR T { iJ ( L ) «-M 2-f ( DD ( L ) /DM ( L ) ) » « 2 ) 

COSPHI ( 1 )=DRR{L)-*ALSD( I )/(!■( 1 ) »D1UL ) ) 

S'l NPH 1(1) =SORT ( 1 . -COSPH 1(1) ) 

ANG = THFTSV+H ( L ) 

(5NFI^(f.OS(AN(j)^Hl(L)+5rN(ANG)«DD(L)/DH(L) ) /F ( I ) * ( - 1 . 1 
TWOR= { COS( ANG )*DD (L ) /DM { L ) “S IN ( ANG ) ^*-0 ( L ) ) /E ( I) « ( -1 , ) 

D\;DSGR ( I ) =ABSD ( I ) * ( rWOR*DwDZSR ( I ) -ONE I DWDZS 1(1)) 

5 DWDSGI ( I )=ABSD( I )^ ( TWOR*DWDZSI ( I ) +ONE I *DWDZSR ( I ) ) 

RETURN 

FND 
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SUUROIITIHE COJGDW (COWZVR, CDWZVI) 

THIG hliBROUTINF CALCULAieS THE COMPLEX CONJUGATE OF THE 
VL'I-OrirV FIELD AT THI ISOLATED VORTFX ‘ 

COMMON / <4 /ALPHA , At* DELTA. DLLTAF » SPANF. XG^F , TANAF , GINA]. SIMA 
12, COSA2. COGA3, COGBT 

COMMON/'S/P] ^ TIa'OPI, PI02, TH<^102» RAD, N, N2, N2X1 , N2P1 
COMMON/6/GVORT , GVORTO, 65(30), 650(30), RS(5]). THtTA5(30), Y5(An 
1), 25(50), D(51), H(52)» RC50), FT0(50), THETA(501, Y{50), Z(‘0), R 
2SV, THET5V, YSTARV, Z5TARV, YV, YVO , 2.V, Z VO , RV, RVO, THETAs' 
COM^'ON/10/LINC 

COMMON/27/DH( 50) » HD(25), DD(50),‘ DRR(50), HH(25) 

CALL DERIVE (YGTARV, ZSTARV, YV,"2V, DVR, DVI ) 
call source (YV, ZV, DVR, DVI, 5VR , SVI ) 

CALL EP5LN (YSTARV, ZSTARV, YV, ZV, DVR, DVl, EVR, EVIJ 
IF (LINC .EO. 0) 1, 2 
I DIDZSR=0. 

DIDZSI = (T. 

GO TO 3 

? call DUINC (YSTARV, ZSTARV, DIDZSR, DIDZSI) 

3 SUMR=0. 

SUMI=0. 

DO 4 K=l, N 
K2 = 2^*K 

YSVMYS=YSTAR V-YS (K2 } 

'^.YSV>^YS-YSTARV+YS(K2 ) 

. /4zSVViZS=Z3TARV-Z5 ( K? ) 

5 YDFNSV= ( YSVMYS-«-*2+ZSVM7.S**? ) «• { YSVP YS**2+Z5VMZS«*2 ) 

E:SUMR=SUMR-GS (K2 )*HH(K) *YS(K2 )*ZSVMZS* ( YSVMYS+YSVPYS) / (TWOPI^YDEN 
^’'SV ) 

4 ;v«5UM I =SUH 1 +GS T K2 ) «HH C K ) -«^YS ( K2 ) * ( YSVMYS^f YSVP Y5-Z SVMZ S**2 ) / ( TWOP I *Y 
-$DENSV) 

|;_SUM R = SUMR+S VR +GVOR'T * E V I / T WO P I +D I DZ S R 

^^;SUMI= - { TANAF-GVORT/ ( 2 . *TWOP I -it YSTAR V ) +SUM I +6V0RT-K E VR /TWOP I -D I DZS I -5 

4 ^ 1 )' 

■ 4ABDV=DVR-i^-»2+DVl*-»-2 
'DZSZVR=DVR/ABDV 
0Z5ZVI=-DVI /A8DV ' 

CDWZVR=DZSZVR*SU[«'tR-DZSZV I^SUMI 
CDWZ V 1 =- { DZSZV I ^SUMR+DZSZVR*SUM I ) 

RETURN 

END- 
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■ SUBROUTINt' (lA^ERV {GVORO} 


IHIS '^iJtiROUT INL' CALCl!LArt -5 THE N!iV.' VORTEX STRENGTH OF THE iSOLATm 
VORTEX BY SA7 ISFy1N(i THF KUTTA CONDITION 

COMMON/4 /ALPHA, AF, DELTA, DLLTAF, SPANP", XSPF , TANAF S 1 NAl , SINA 
12, COSA2, COSA3, COSBT 

COMMO.N/S/PI , TWOPI, PI02, THPIO?, RAD, N, N?» N2M1, N2P1 
COMMON/G/GVOKT, GVORTO, GS(SOi, GSO(FnI, R5i<3\), THETAS(Sf>), YS«‘^0 
]), ZS(SO), D(S1)» H(52), K(50), R0(50), 1HI-:TA(50), YISOI* 2(30), R 
2SV, THFTSV^ YSTARV, ZSIARV, YV , YVn , ZV , ZVOf RV, RVO, THl:TA\!- 
COMMON/IO/LINC 


SUM = (^. 

DO ] 1=^3, N 

I2 = ?^f I 

1 5UM = SU-M+ (H( I2 + ?) -H( 12-2 ) )*GS( T2 ) «COS ( THETAS ( 12 ) ) /RS ( I? ) 

IF (lINC .LO. 0) 2, 3 

2 D1DZSI=0, 

GO TO A 

3 CALL DWINC (0., 0., DIDZSR, DIDZSI) 

A DIV=COS(THETSV)/(PI»RSV} 

GV0RD = TANAF-GS(2 ) *H ( A ) **2*C0S C THFT AS ( 2 ) ) / ( TWOP I * ( H ( A ) -H ( 2 ) )«RS(2 ) ) 
1-GS{A}*(H(6)-H(2)*H(A)/{H{A)“H(2 ) ) ) «C0S { THETAS ( A ) ) / { TWOP 1 5fRS ( A ) ) -S 
2UM/TW0PI+DIDZSI 
GVORO=GVORD/DI V 
RETURN 
END 
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sun.K’oun^iF lovOPI (istopd 

IHlh SUBROUTINE IS UFSIGNEU TO SATISFY THF PRISSURE CONDI nOM 
CQ'MNOM/?/DX , XC« X, XXO , 

COMMON/5/PI? 1W0C>1. PJ02»- THPl’02» RAD» N» N2» N2M1<. N?d 1 
.-COMWON/S/GVORT , OYORTO, GS(50), GSG(Sn), RS(5l-)» THFTASISrn. 

:T ) » ZS(SO)i 0(A]), H(S2 )vR(50)» RO(SO), THETA { ),' Y ( SO ) , R 

• 2SV, THFTSV» YSTARV, ZSTARV» YV, YVn , ZV » ZVO» RV-. RVO» THFlAV 

COMMON/ll/DWDSGR(2S) » OUDSGI(2?)» E(2S)» COSPH I ( ? 5 ) , SINPHU2S), A 
1BSD(?5) 

.COMMON / 18 /ACL- 

COMMON/24/NLOOP1 , ACC1» SUM1(2S), SUM10{2S)» SUMIUZS), SSUMK2S). 
1 SUM2(25)» SUM20(2S), SUMZKZSJf SSUM2(2S)» CSPIXO(2S)‘, 'D‘vOSXO(2S) 
2» RX0(2A) 

COM'^ON/27/DH { 50 ) > HD (25), DD(50), DRR(50)* HH(?S) 

DIMENSION AD (25), 6(50), ER(25) 

I STOP 1=0 
ACC=ACC1 
NITFR=C 

CALL COFTGR ( 1 ) 

6 IPASS=0 

DO 3 1=1 , N2 

3 G(I)=GS(I) 

GV0RD=GV0RT 
CALL USEFUL 
DO 1 1=1, N 
L=2«I-1 ’ 

SUMll ( I ) =CnSPHI ( I ) 

SSUMK n=SUMl ( I ) + (SUMlO( I )+SUMll { I) ) « ( R ( L ) -ROC L ) ) /? . 

SUM21 ( I )=DWDSGR( I -) 

SSUM2C I )=SUK2 ( I ) + (SUM20( I )+SUM21 ( n }* (X“XO)/2. 

AD( 1 ) = (CSPIX0( I )*RX0( I l+SSUMl-DWDSXOl I •) »XX 0-SSUM2 ( I ) )*ABSD( I )/F.( I ) 
SUM=SSUM=0. 

DO 2 1=1, N 
IN = N + 1-I. 

SUM=SUM+SSUM 
L = ?-»IN-1 

GS{L)=(-GVORT-SUM)/( .5*DH(L)+AD( IN) ) 

2 SSUM=HD(IN)*GS(L) 

CALL C0FI6R (0) 

DO 4 1=1, N 
I2=2«I 

ER { I ) =ABS ( ( GS ( 12 ) “G ( I 2 ) ) /G ( I 2 ) ) 

• IF (ERd) .GT. ACCl) IPASS=1 

4 CONTINUE 

- CALL gamer V (GVORT) 

ERROR = AB.S( ( GVORT -GVORD ) /GVORD) 

IF (ERROR .LL. AlC.) .AND. 1 PASS .EO. 0) GO TO 5 
MITER = N1 TER+ 1 

IF (NITI.R .LL . NL00P1 ) GO TO o 

ACC 1 =ACC 1 +ACC 

IF (ACCl .GT. ACL ' ) 22, 

22 IST0P1=1 

PRINT 20, fFP(I), 1=1, N), ERROR, ACCl 
RETURN 

R GVORT=(GVORT+GVORD) /2. 

DO 10 1=1, N2 
in GS( I ) = (GS( I )+G( I ) )/2. 

GO TO & 


41 



fNi Tv: 


COM T I Ml !!^ 

ir (Acn .Ko.-- ACC) 00 T.o 

PRINT 20, ( FR.{ I } » -1 = 1 , N), ERROR, ACC] 

ACC]=^ACC 
A RETURN 

a FORMAT 2uHERR0Rb ARE TOO LARGE, 2X, 3HE R = 6 ( El 2 . '5 , 2X )//[;<, 

inERROR^ E]2.I>, 2X, 3fiACC I=E 12 . ) 

END - ' ' • - 
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SUBWOUTlwr LOC;P2 (iM» ISTOP?) 


Ttii'’) runRouTii'jr: is ijrsiGNFi) to satksfy the: force; halancf i:ouatfom 

XEi, X. XXO 

C0'1MC)N/5/PI » T«'0pr, no?, -IHPIC?, pad, N, N?, H2M3, N2P] 
COMMOH/fi/CiVORT , OV/GRTO, GS(50)* GSOEbO), RSE'"'!), THET As ( fl ) , YE('>0 
1), 2S(50), D(S1), H(E>2), R(50), R0(5C), THETAEAO), Y(50), Z{'=>0), R 
2SV, THETSV? YOrARV, ZSTARV, YV, YVO» ZV, ZVO, RV, RVn , THcTAV 
COMMON/ 18 /ACL 

COf1MON/25/NLOOP2 , ACC2» STEP 

Dlf/ENSrON Y5V(S), ?SV(S), V0RLFT(4)» ZFTAR(A), ZfTAT<Al, FR{4.) 


!5T0P? = t) 

ACC=ACr? 

IF (H) 3, ?., 3 
? CALL LOOP] (ISTOPl) 

IF (I STOP 1 .FQ, 1) GO TO ?8 
RFTDPN 

in PRINT 31 , ER ( AI , ACC? 

?8 ISTnP?=] 

RETURN 
3 NITER=0 
1=1 

YSV( 1 ) =Y5TARV 
2SV ( 1 ) =2STARV 
27 CALL LOOPl (ISTOPII 

IF (JSTOPl .EO. l‘) GO TO 10 

NITER=NrTFR+l 

VOPLFT ( I 1 =GV0RT-H-YSTARV*4. 

call COJGDW (CDWZVR, CD'VZVn 

ZETAR( I ) = (GVOI^T-GVORTO) /DX*(YV-Y(N2 ) )-GVORT« ( CDWZVR- ( YV-YVn } /QX ) 
ZETA] ( I) = (GVORT-GVORTO 1 /DX* { ZV-Z ( N2 ) 1 -G VORT* < CDWZ V I - ( ZV-ZVO ) /DX ) 
ERE I inZETAR ( I )«-*2 + ZETAI E I ]*^i-2) /VORLFTE I ) 

IF EFPEl) .LE. ACC2) 1, A 
TF (ACC? .FO., ACC) GO TO 35 
PRINT Al, ERE4], ACC2 
ACC2=ACC 
35 RETURN 

A GO YO ( 5 , 6 , 7 , 8 ) I 

5 YSTARV=YSVE 1 ) +STEP 
YSVE?)=YSTARV 
ZSVE2) =ZSTARV 

I=? 

GO TO R 

6 YSTARV=YSVE 1 1 
Y5VC3)=YSTARV 
25TARV=ZSVE1 )+STEP 
ZSV(3)=ZSTARV , 

1=3 

GO TO 9 
9 PSVOLn=PSV 

RSV = SQRT ( YSTARV«^«-2+ZSTARV<f^^? ) 

DLLPSV=PSV-RSVOLD 
TERM = DELRSV/HEN2 ) 

DO ?6 J=l» N 

?6 DEI ? ) =DE I ? )+DFLRSV-TERE^'«HE 1 ? ) 

GO TO 77 

7 7fRMl = Zr:TARE 2 ) *Z FT A I ( 3 ) -2 £T AR C 3 ) *-Z£ TA I ( ? ) 

TFR^V=ZrTARE 1 )*ZFTAI ( 3 1 -ZET AR ( 3 ) *ZETAI ( 1 1 
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1 IRM'^t^Zr-TAP ( 1 ) rfZt 1 AR ( ? ) -Zf-TAR{ 7)-)fZETAn 1 ) 

DLTFR" -TTRvi-rFRMA 

YbV/ = { YbVi .1 ] *n rR!'a-YbVC ? j »-TtRM2+Y<>V( 3 )-»T{;RM3) /DF.TEKM 
ZS\/Z=(ZSV (1 ) * reR-U-ZSVt ? l^i-TERM?+ZSV( 3 )*Tf:RK31 /iJt rf-KM 
IP IAf3,^( ( YSV2-Y6TAPV)/r5,rARV) .GF. .5) 12. 13 

13 IF IAB5( (ZSVZ'ZG1APV)/7.STARV) ,GF. .3) 12, 14 
P ^PR"AX=A^^'^>'1 I FR( I ) , LR(?), FR(3)) 

IF ftR(A) .Gl. ERRMAX) 15, lA 

15 IF (MTRY .GT. 6) 30, IR 

30 ACC2=Af C2 f-ACC 

IF (ACC2 .GT. ACL ) 10, IB 

18 NTRY^MTRY+1 
iNITCR = M] TER-I 

ERP-MIN = AMIN1 CER( 1 ) , £R(2). ERf3)I 
on l<9 J=1 , 3 

IF fER(J) .NF. E’RPMIN) 10, 2 0 

20 L^.J 

GO TO 21 

19 CO^mNUF 

21 VSV2=YSV(L )+ C YSV2-YSV(L ) J/2. 

ZSVZ=ZSV(L)+(Z5VZ-Z5V(L) 1/2. 

14 YSTA.RV=YSVZ 
Z6TARV=ZSVZ 
1=4 

PM = 1 
GO TO 9 

16 IF INITFR .EQ. NL00P2) 3U, 22 

22 NTRY=0 

DO 23 J=1 , 3 

IF (ER(J) .£Q. ERRMAX) 24, 23 

COMTINUF 

YSV( J}=Y5VZ 

ZSV{ J) =zsvz ■ 

ZETARI J)=Z£TAR<4 5 
ZETAI ( J)=2ETAI <4) 

ER( J)=L'R(4) 

VORLFTI J)=V0RLFT(4) 

60 TO 7 

TAU=ATAN( (ZSVZ-ZSTARVI/(YSVZ-Y3TARV) > 

YGVZ=YSTARV+.01*COS ( TAU ) 

ZGV2=ZSTARV+. 01 ^-SIN ( TAU) 

GO TO 14 

31 FORr^AT (IHU, 18HERRUR K5 TOO LARGf , 2X. 6HERR0R=M 2 . S » 2X, 
1F12.5) 

END 


5HAC C^ 
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r-.UBROl)T I.«)F. L00P3 (ISTOP3) 


C THIS SUHROUTiMf. IS UIESJGH!;D 10 SATISFY THF NORMAL VELOCITY 

C0MWN/1/RADIUS» AZAXIS* BYA/ISj S» OLTA 
C0'^’'<''‘N/3/DX > XT, X, XXO 

COMMON/'3/r-I , T'OPI, P102, THP102, PAD, N, MP, N2M1, N?P1 
COMMON/6 /GVOR T , G'/ORTO, O5!50), GSP{Sn), RS(Sl), TMCTT S < 6 ) , 

II, ZS(50), D(S1), H(62), R ( Sn ) , JTO ( SO I , THETA('O), Y(SO), / ( '=■0' , R 
2SV, THFTSV, YSTARV, ZSTARV, YV, "Vvo, ZV, ZVO. RV, RVO, THETAV 
COMMON/l 1 /DWDSGR (25 ) , UWUSGI(25I, t(25)» COSPHI<75J, SINPl^^(2^), A 
1HSD( 25 ) 

COMMON/ IB/ ACL 

COMKON/26/NLOOP3, ACC3, COTPHUZA), THETA0{50) 

C0^i^'i0N/27/DH( 50) , HD(?S), D0(50), DRR(50), HH ( ? S ) 

DIMENSION FTA(25), D0LD(25), DI:l(?5), ER(25) 

ISTOP3==0 

ACC=ACC3 

delold=o. 

APR I''E=. 05 
N ITFR = 0 
M=0 

1? call loop? (M, I5T0P2) 

IF (1ST0P2 .FO. 1} GO TO 23 
DO 6 1=1, N 
6 D0LD( I )^D( 2* I ) 

DO 1 1 = 1, N 

y l =2 * I - 1 

% FTA( I ) = ( (R(L >-R'HL)-R0( L )*COTPHI ( I )* ( THAJA (L )-TMFTA0I L ) ) ) /DX-PwnSG 
'T-1 n I ) /STNPHI ( I.) ) /DWDSGRI I )*SI MPHI ( I ) 

1(£ ETA( I )=AS!N(ETA( n ) 

] DCLM )=-FTA(1)*(D(2)»*? + RSV*«?-?.«D(2)*RSV*C0S(DH(1 )) ) / ( RSV^SI N ( DH 
1) ) ! 

DO 2 1=2, N 
f L=?«I-1 
i’ I2=?«I 
? I2M?=I2~7 

2 ' DEL ( I ) = (0( I? )*DEL( T-1 )1 /DC I?M?) -ETM I)»{D(I?)^fK2 + 0< I?M2 ) I 

1? )*D( I2M?)«C0S(DH(L) ) ) / ( D ( I 2M? ) ^^SIM ( DH ( U ) ) 

DO 3 1=1, N 

3 EP( n=AB5(DEL( ! )/D(2*I ) ) 

N ITFR = WI TF.K + 1 
ERRMAX=ARRAYMAX( ER, N) 

IF (FRRMAX .LE. ACC3 ) 7,8 
8 IF (NITFP .GF. NLOOP3) 9, 10 
P ACC3=ACC3+ACC 

ir (ACC3 .or. ACL ) 13, 10 

13 PRINT 2S. LRRMAX, ACC3 
?3 isroP3=i 

Rf torn 
10 M = 1 

IF (DEL (N)-»DFLOL0) 4, 5, 5 
5 DElOLD=DFL (M l 
GO TM lA 

4 APRT'''E = AdrIwF/7. 

DELOLD=0. 

14 IF (F.RRNAX ,LE. APRJME) 15, 16 

15 FACTER=1, 

GO TO 17 

16 FACTER=APRiME/fcRRMAX 
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17 fONTPI--'- 

00 m 'J 

IB 0( n -DnL0( I l+OfLt I ) 

G(' TO IE 

7 If (ACTB '.CO^ ACC) GO TO 00 
PRIMP ?'j), ERRMAX» ACC3 
ACr^ = ACC 

30 RETURN 

?6 roi^MAT (IHO, IRHERROR IS TOO LARGE, 2X, 6HEPNAX=El ? . 3 , ?X, 
1C12.5) 

END 


SUACC'’ 
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■^>UBRniiT[ME LfFT (1C) 


THIS SntJRGlJTlNi: CA| Cl'LATfS LINFAR LIFT CO£F. AMO ''ASONL I ME AP 1. f( T 
COEFF.. < = f, CALCULATES EVEf^YlhlFG AT XO STATION. 

K=l» CALCULATLS EVrRYTHING AT ANY STATION. 

COMMON/i/RADlU^, A21AXIS, BYAXIS, S, BETA 
COAMDN/-H/DX, X> . X, XXO 

COMMON/ A/ALPHA, AL, DELT/' , *DF LTAF , SPANF, XSPF , TANAF. SINAI, .^INA 
12, COSA2, C0SA3, COSBT 

COMMOM/E/PI, TlOPI. PJ02, THPI02, RAD, N, N2, M2M] , N?PT 
COMMON/6/GVORT » GVOFUO, GS(50), GSFHFD), RSCSl), THFTAS(SO). YSi'^r 
1). ZS(S'T), D(F1), H(5?), RJ50), RO('^O), THETA(AC), Y(SO), Z('-O), R 
2SV, THnSV, Y^TARV, ZSTARV, YV , YVO,‘ ZV , ZVO. RV, RVO, THETAV 
CCMMON/8/SPRSQS, SMRSOS 
COMMON/S/FPS 
COMMON/] 0/LI NC 

COMMON/27/D1-H S(i) , HD{?5), DD(AO), DRR(Sn), Hh ( 2 A ) 

COMMON/Pn/SUM? , SUM'^0, SUM31« SS1JM3, SUM4 , SUMAP, SUM41 , SSUf'^4 , c;t. 
IMS, SUM5F-, SUM51, SSUM5, SUMS, SUM6 0, SUM61, SSUM6 
COMMON/30/TLIFT, UPi^FR, XLOWFR 

IF (K .FO. 0) 5, 6 

5 SUM30=S 
SU‘'*3 = .5>‘S*X 

CL1=PI*SINA1*C0SA2^«-TLI FT/ (C0SRT*StJM3 1 
SU'>' = 0, 

DO 40 1=1, N 
L = 2 

40 SUM = SUM+GS(L >*2.*YS(L)*HH( I ) 

CL2 = 2.* tGVORT»2.-K-YSTARV+SUM/2. ) *C0SA3 / ( COSBT^fSUM3*2 . ) 

CL=CL1+CL? 

SUM4D=CL1*2.*-SUM3 
SUMSO=f L2*?.*SUM3 
SUM4=2 .*CLl/3.*2.*SUM3*xn 
SUM5 = 2 .*CL?/3.’^2 .*SUM3*X0 
AOM] =SUM4/(2 .-»sSUM3*X0 ) 

AOM2=SUM5/(2 .*SUM3*X0 ) 

A0M=A0MI+A0M2 

RATIOO=AOM/CL 

CLIN'O. 

SIJM60=0. 

SUM6=0. 

SSUM6=0. 

GO TO 8 

6 SUM31=S 

SSUM3=SUM3+(SUM30+SUM31 )*(X~X0)/2. 

CLl = PIfrSINAl*C0SA2<^TL I FT / ( C0SBT*SSUM3 ) 

5UM=0. 

no 1 I=I , N 

L = 2»I 

1 SUM=SUM+nS(L)*2.*Y5(L)^HH( I I 

CL2=?.*(GVORT«2. «YSTARV+MJM/2. ) *C0SA3 / { COSBT«SSUM3«2 . ) 

Cl =CLl < C L2 

SUM41 = CLI «2« «SSUM3 

SUMS1=CL2»?. «SSUM3 

S5um4 = s(;M 4+( x+xo ) X ( SUM4 1 -sn'-’An ) n . 

A OM 1 -S SUM4 / { ? . « S SUM 3 *•>' » 

S5UMS=SUM''+ ( X+XO ) X ( SUMSl -SUMSO } /? , 

A0M? = SSUM'^/( 2.XSSUM3X-X ) 

A0M=A0M1+A0M? 

RAT 100=A0M/CL 
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IF (LINT ,E0. 0) 2 

^ CMN^O. 

GO TO 7 
? R2 = PAD 

THFT=-ATAN2<S«-^f2-R2» 2»*RADIUS«-R ) 

IF (THFT .LF. 0.) THFT==THFT+TWOP I 

CL FN-FPS*COSA^i«'( ( TL I FT+R2 ) *AC05 ( 7 . -»R AD I US^ i>/ ( S«« 2 + R2 ) ) +-7 4^7^^: SPR .0 
7S«!-2-7LlFT-R2 ) +SMRSOS *«■ ( P I ’^i.MRSOS/2 . -4 . *RAD lUS ) )/ ( COSB 7*SSUM3 ) *7 . 

CL=CL+CLIN 
SUM61 =CLIN'*2.*SSUM3 
550*^6 = 5006+ ( X+XO )*{ SOM61-SUM60} / 7, . 

A0‘^T^SSuV6/{ 2.+S5UM3 + X) 

AO'1=AOM+AO''1 I 

7 S0M^0=RUM3I 
. 5UM40 = SO.'141 

SUM50=SOM51 
SUM3=5SUM3 
SUM4=SSUM4 
SUM5=5S0M5 
SUM6 = S5U'>^6 

8 PRINT 4, CL1» CL7, CLIN, CL, AOMl, A0N2 , AQM, RATIOO 

4 FORMAT (IHO, 12X , 18HLINEAR LIFT COEFF=E12.5, 2X, 21HNONLIKFAK L IF 
IT COEFF=F12.5, 2X, 17HINC0T LIFT C0EFF=E1? . 5 //33X , 17HT0TAL LIFT c 
20EFF = E12.5//13X, 14HLINEAR MOMENT=E 12 .5 , 2X» 1 7HN0NL I NEAR MCOENT-E 
312.6, 2X, 19HT0TAL MOMENT COEFF^E 12 . 5// 33X , 6HCM/CL=F 1 0. 5 ) 

RETURN 

END 
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M-fVxOMT 1 '1 P'?f IT'S I L'M 


’lil£. '.V 1 P'H 1 1 '4L PiaPAR' 5 CALaM-ATLS Thf; r> r ^;^.pr , 

PRt: 3 Ai-L TlOfir, M nr Ffi'31 c-TAriC'U 

(. yL-ruL \TF3 rnr. PR^3oU^^. at st coni^ siat:o-\. 

C Ok/ : /RADlIJi-.. a; 1 3 •. r*VA<li>» St BPTA 

f M /'’/ MSHAPI 

(; f, IMV-Jti/T/Ox. >XXt XXO 

c., V/04/A/PI . liOPi, PI02, TnPHV* RAft R- N ’ - N?Mj , m 2P 1 
COM' 0 4/6/GVORT j GVCRTO, CS(5Plt GSO{50)- RS(kl), Thf- T'' S { 30 ) t 
1), ZSC'Mt 0(31), H(32), R(3r!, -0(30), TH'-''A[kO), Y(‘‘-0), '(kO), P 
?3V, UriSV, YriARV, ZS7ARV, YV, AVO, ZV, ZVO, PV, PVO. THCTkV 
Cn<'MOM/i''/x, MPRf, inojT 
CO’ '^0.M/^''/TLIFT, CPPFR, XLOOCP 
C Jf' 'OM / 3 j /CZV;> , CZ k ( 2 3 ) 

YPS(2v)2), ZPS(202), YfMP'^?), ZP(202), PHIX0(202)-. PH 1 X < 2 
1U2), PRk(202)' NOPHI(202), PHIXY(2C2), PHIXZI2Q2) 

COMOlEX CZV3, CZS 

N<lr:MpRr^] 

,M2 = N!PR{>2 
M3 = 2»,MDPf4-? 

DO 3 KJ=1, M3 
V. NODhI (< J)=0 
■'.MF (LM .EO» 1) 1, 2 
i;‘-' X=X+DX 
XXX =X 

CALL NCONSH (3 ) 

S CALL SETUP 
YPS n )-0. 

YZ YPS (Ml ) = 0. 

'’>= YPS(M2)=0. 
fs YPS(M3)^r, 

/ ZPS ( 1)=.OOU3*UPPEP 
’* ZPS (Ml ) = , 99097i-UPPER 
ZPS ( M? ) - , 0001 *XLOWF.R 
?PS [M3)=.9999^XL0WLR 

CALL YZCOMP (YPS( ]), Z^St 1), YP ( 1), ZP ( 1)) 

CALL YZCOKP (rPS(Ml), ZPS(Ml). YP(MD, ZPlMD) 

CALL YZCOMP (YPS(i'!2>, ZPS (M2), YP(M?), ZP(M?)) 
call YZO ip (YPS(M3), ZPS(M31, YP(M3), ZP(M3)) 

UPPER-UPPER//. PPE 
XLOWEP = Xt OWER/OPRE 
DO 44 1 = 2, A!PRE 
YPS ( I )=n. 

/PS [ I ) =UPPEP* ( I “1 ) 

CALL YZCOMP OPSII), ZPS(I), YD(I), ZP ( I ) ) 

J= 1 +M1 
Y(^S ( J) = 0. 

ZPS ( J)=Xl OWFR^-( I-l ) 

44 CALL YZCOMP (YPS(J), ZPS(J), YP(J), ZP(J>) 

X=X-DX 
XXX=X 

CALL NCOOSH U) 

CALL SETUP 

CZVS=CrPL X (YSTARV, ZSTARV) 

DO ST I=], N 
I2 = ?k I 

SI CZS ( 1) =''^1PLX ( YS( 12 ) , ZS(I?)) 
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DO t = l* Ml 

call Y7'-,COf> (YP(I)» /P([)5 Y»S(I), ZPb(I}) 

IF (7PS(I) .rc. 0.) /Pi,( I )=ZPS( I )+l.!--8 

CALL Phil (vp.(I). 7Pi>(n, YPtI)» ZPU)» PHIXO(I), ML'’) 
ir ML'- .LO. I) AS, SL 
AS MOPHlin^l 
J=I+M1 

CALL YZSChp (YPfJ), ZP{J), YPS(J), ?PSfJ)) 

IF (((AZAXIS 0) .AND, (NSHAPE ,FQ. 1)) .OR. ((RADIUS .FO. v'' > 

l.AND. (MSHAPL .E(^. 0))) 59» 60 
SO 2PS ( J) =-7PS( J ) 

GO TO 40 

60 IF (ZP(J) .EL. 0.) 2P5( J)=-ZPS( J) 

IF (2PS(J) .EO.0.1 2PSS J)=ZPS( Jl-l.C-B 
AC) CALL PHir (YPS(J), ZPS(J). YP[JL. ZP(J), PHlXn(J), MLP)t^ 

IF (MLP .FQ. 1) 61. 46 

61 NOPNI(J)=l 
46 CONT I NUF 

RETURN 
2 PRINT 61 
YPS < 1 ) =-{). 

YPS (hi ) =o. 

YPS(MZ)=0. 

YPS CM3 ) =0, 

ZP5 ( 1 ) =,0001*UDPCP 
ZPS (hi ) = .9999^<UPPER 
ZPS(M2) = .0001*XLOi-'ER 
ZPS (M3)=.9999»XL0WER 

CALL Y2C0MP (YPS( 3), ZPS( 1), YP ( 1), ZP ( D) 

CALL YZCOMP (YPS(Ml), ZPS(Ml), YP(Ml), ZP(Ml)) 

CALL YZCOMP (YPS(M?), 2PS(M2), YP(M2>» ZP(M?)) 

CALL YZCOMP (YPS('lo}, ZPS(M3). YP(M3)» ZP(M3)1 

UPPER=UPPER/NPRE 

XLOWER^X LOWER /NPRE 

DO 50 I-?, MPRE 

YPS ( I }=0. 

ZPS ( I 1 :^UPPER« ( 1“1 ) 

CALL Yzcni^p (YPS(n, zPS(i). YP(i), zp(in 

- J=I+N'l 
YPS ( j) =n. 

ZPS (J) =XLOVFR<^ t I-l ) 

SO CALI YZCOMP (YPS(J). ZPS(J), YP(J), ZP{J)) 

CZVS=CMPL<( YSTARV, ZSTARV) 

DO 52 1=1, N 
I2 = 2« I 

52 CZS ( I] =CMPLX (YS( 12 ) , ZS(I2)) 

DO 53 1=2, Ml 
J=I+M1 
IM1=I“1 
JM1=J“1 

CALL PHI I (YP5(n, ZPS(I), YP ( T) , ZP ( D , PHIX (I), MLP ) 

IF (,MLP .EO. I) 54, 55 
5 A NOPHI ( n = l 

55 call phi I(YPS( IMl ) , ZPS I I) , YP(JMl), ZP(I), PHIXY(I), MLP) 

IF (MLP ,E0. 1) A, 5 

A NOPHI ( 1 ) = 1 

5 CALL PHTr(YPS(I), ZPSIIMl), YP ( I ) , ZP(TMl), PHIXZd), MU’) 

IF (MLP ,EQ. 1) 6, 7 

6 NOPH n n -- 1 

7 CALL PRF'.S (YI’S(I), ZPS(I), YP ( I ) , ZP ( I) , PHIXO(IMl), PHlX(I), PHI 
1XY( 1 ) , Pm IXZ ( I ) , PRE ( I) ) 

CALL Pflll (YP'.(J), ZPS(J), YP(J), ZP(JI, PHIX (J), MLP) 
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' ir (*<LP »ro. 1 ) ‘56, hi 

'■'6 noPHTfJl^l 

67 CALI. PMl I ( YPS: JM.l ) • ZPSM). YP{Jf'l), ZP { J ) , Pl-i!XY(J). 

IF fWLP .‘“Q. 1) a» «5 
P /‘^0P!IHF1 = 1 

9 CALL PHII(YP6(J), ZPS(JMl), YP(J), ZPfJMl), PHIXZ(J), MLP ) 
jr (MLP .ro, 1) ivO, 11 
] n NOPin ( i ) -1 

n CALL PRESS (YPS(Jt, ZPS(J), YP(J). ZP(J)» PHIXOtJMn, PHIX(J) 
IXYf J) , PHIXZ ( J1 , PPE( J) ) 

PRINT 6P, I, YP(I). ZP(I), PRS(I), NOPMJ(I), T, YP ( J t , 7P(J), 
IJ) , NOPHI ( J) 

DO 12 KJ=1, M3 
12 NOPHI (KJ)=0 

IPRTT=IPR [T+1 

61 FORMAT (IHu, IX, 36HPRESSURE COCFF. ON UPPER SUREACC ART, 24X 
IPRESSURE COEFE. ON LOWE» SURFACE APE //IHO, IX , IHI, 8X, IHY, 
2 1H2, 12X, 2HCP, 7X, IH*, IX, IHS , IX, IHI, 8X, IHY, 13X, IHZ 
3, 2HCP, 7X, lH«-/flX. 1H£) 

60 FORMAT (SIX, IHS/IX, 12, 3(2X, E12.5), 2X, II, IX, 1H$, 1X» I 
12X, F12.5), 2X, ID 
RETURN 
END 


Phil 

PRF( 

36H 
1 3X , 
12a 
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‘^UBRCIITINE PHIKYPS. /PS^ YP, ZP» PHI, M) 

THir. ^IJbROUTIfJE CALCULATES ThF IMAGINE PART OF THE POTENTIAL 
font I ON AT ANY P05 I Tl ON X 

COWPON/] /RADIUS, A2AX1S, BYaXIS, S, BETA 
C Of MO ^ I / 2 / N S H A P f: 

COMMON/A/ALPIIA , AE, DFLTA, DFLTAF, SPANF, XSPF , TANAF, SINAI, SIN^* 
12, COSA2, COSAJ, COSBT 

COMMON/5/PI, TWOPl, P102, THP]02, RAD, N, N2, M2M1, N2P1 
C0MM0N/6/0VORT, GVORTO, GS ( 50 ! , 'GSO ( 50 ) , PS(51)» ThtTASIfjO), YSIAf', 
]l, ZS(5'.), DC 51), HC52), RC50), ROC 50), THETA! 50 I, YCSO), Z C 5v'i ) , R 
2SV, THFTSV, YSTARV, ZSTARV- YV. YVO , ZV, ZVO« RV, RVO, THETAw 
COMMON/15 /ROOTR » ROOT I 
COMMON / 1 P/DBDX , DRODX 

COMMON /27/DHC 5 )) , HDC25), 0DC5O), ORRC50), HHC25) 

C0MM0N/31/CZVS, CZSC25) 

COMPLEX CZPS, CZVS, CZS , CL 

M=0 

CZP3=CMPLX( YPS, ZPS) 

SUM=0, 

DO 1 1=1, N 
12 = 2* I 

CL = CLOGC CCZP5-CZSC I ) ) / C CZPS+CON JG C CZS C I ) ) ) ) 

IF CZPS .GT. ZSC 12 ) ) 2, 3 

3 IF (ZPS .EO. ZSC 12) ) 4, 5 

5 IF (ZPS .GT. 0.) 6, 7 

6 .|l=CL+( TWOP l* ( 0. » 1.}) 

:go to 2 

7 J:F (ZPS .EO. 0.) 4, 2 
2 f'ACT=AIMAG(CL) 

1 SUM=SUM1+GS { I 2 )*HHC I )*FACT 

OL=CL06( CCZPS-CZVS) / (CZPS+CONJGCCZVS) ) ) 

>i;F CZPS .GT. ZSTARV) 8, 9 
9 *’^F (ZPS .EQ. ZSTARV) 4, 10 

10 fjf CZPS .GT, 0.) 11, 12 

11 ,eL=CL+C TWOPl* CO. , 1.)) 

Go TO 8 

12 IF (ZPS .FO. 0.) 4, 8 

8 FACT = AIMAG (CL ) 

IF (NSHAPE .FQ. 0) 1a, i5 

14 WSR=RADIUS*DRDDX*.5*ALCG( YP**2+ZP**2 ) 

GO TO 16 

15 WSR=A'ZAXI5*DBDX *. 5*AL0G ( ( YP^ ROOTR ) **2+ ( 2 P+ROOT I ) *“2 1 

16 PHI =TAMAr*ZPS+GVOPT*FACT/TWOPI + SUM/ (2.«TW0PJ )+V'SR 
RETURN 

4 M=1 
PHI=0. 

RETURN 

END 
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sufjpoi^n iNP pRc:>s (yps, zpp, yp, phjo» phi » phixy. phik?, cp5 

C THIS 51i!»P0UTlML CAl CULATf f> PRrSSlJPE COEFFICIENT 

CO.MK'Or-i/3/nx* XO, X, XXO 

common /A /AL f'HA » At, DELTA, DLLTAF, SPANF, XSPF, TAj'JAF. SINAI, SINA 
12, C05A?, COSA3, COSBT 

COPMON/3/P1, T-OOI. PI02, THPIO?, RAD, N, N2, NFM], N2P1 
COKMON/A/AVORT , GVORTO, CM 5 0 ) , GSO(^O), RS(S]), THCTAAl^o), VSfAO 
1), ZS(A'’i), D(S1), HfS?). R(S~i), RO(SO), THFTA(Sn>, Y(AO), R 

2SV, THFTSV, Y'^TARV, 2STARV, YV,-*YVn, 2V, 2V0, RV, PVO, T'lETAV 
COMMON/27/DH (50) , 1.0(25), DD(50), DRR(50), HH(25) 

CALL DERIVE (VPS, 2PS, fP , ZP, DR, DI) 

ABSDS0=1. /(DR*fr2+DI ^^^2 ) 

CALL SOURCE (YP, 2P, DR, DI , SR, SI) 

YSMYSV^YPS-YSTARV 
YSPYSV=YPS+YS7 APV . 

2SMZSV = 7PS-ZS1 ARV 

YDENV- ( YSMYSVir«? + ZSMZSV*-«2 ) » ( YSPySV«« 2 + 2SMZSV*^2 > 

SUMP-0. 

SUM 1=0. 

DO 2 K=l, N 
K2-?«K 

YSMyS-yt>S-YS(K?) 

YSPYS=YPS+YS(X2) 

2SMZS=ZPS-ZS(K?) 

YDENS= ( YSMYS»*2+ZSMZS*»? )^ ( YSPYS**2+ZSMZS**2) 

SUMR=SU.MR+6S ( K2 ) *HH( K ) *YS ( K.2 ) 55-ZSMZS* ( YSMYS+YSP YS ) / { Tvv'OP I^YDFNS ) 

2 SUMI =SU.MI +GS ( K2 ) *HH( K ) *YS ( IA2 ) * ( YSMYS*YSPYS-Z5MZS**2 1 / ( TWOP I -^YDr NS ) 
DWDZSR =-GVORT *YSTARV*ZSMZSV ^ ( YSMYSV+Y5P YSV ) / ( P I * YDENV ) -SiJ.MR + 
ISR 

DWDZSr =-(GVORT *YSTARV* ( YSMYSV*YSPYSV-ZSM7 SV-k-*2 ) / ( P I »YDFNV ) +S 

lUf^I+TANAF-SI ) 

Dl^SQ^D''DZSR^s^f-?+OWDZSI «-*2 
PHIX=(PH1 KY+PrilXZ-PHI-PHlfi) / ( X-XO ) 

C f? = S I NA2- ( 2 . «PH r X+ABSDSQ*DWSQ ) *COS A2 ' 

RETURN 

END 
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